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1. UVOD 


Da bi se elektroničko digitalno računalo prilagodilo za proračunavanje dinamičkih susta- 
va i da bi se način simuliranja procesa i programiranje približili odgovarajućim postupcima kod 
rada s analognim elektroničkim računalima, razvijeni su posebni programi za digitalno računalo. 
Jedan od takvih programa kojeg ima digitalno računalo Elektrotehničkog fakulteta Sveučilišta u 
Zagrebu (IBM 1130) je 1130 CSMP (CONTINUOUS SYSTEM MODELING PROGRAM). 


Ovaj program vrlo je prikladan i praktičan (i jednostavan), jer nije potrebno poznavati 
niti jedan programski jezik digitalnih računala, a programiranje je vrlo jednostavno i za one ko- 
ji nisu i ne žele biti programeri. 


Programiranje je u principu vrlo slično programiranju za analogno računalo. Blokovsku 
shemu zadanog nekog, na primjer, regulacijskog kruga potrebno je nacrtati pomoću elemenata koji 
stoje na raspolaganju kod programa 1130 CSMP. Pored 25 raznih standardnih elemenata postoje i 
specijalni elementi koje se može definirati prema posebnim potrebama. Maksimalni broj elemenata 
koji se može upotrijebiti kod rješavanja jednog problema je 75. Od toga može biti najviše 25 in- 
tegratora, ali mora biti upotrebljen barem jedan integrator, jer je to u stvari program za rje- 
šavanje općih diferencijalnih jednadžbi. Jedna od prednosti pred analognim računalom je i ta da 
nije potrebno tzv. skaliranje varijabli i da se mogu dobiti numerički podaci iz bilo kojeg mje- 
sta blokovske sheme. 


Svaki element prikazan je svojim simbolom u obliku grafičkog znaka i sa slovom, odnosno 





cznakom. 


Nakon što je pomoću raspoloživih elemenata nacrtana blokovska shema sa svim potrebnim 
veličinama i vrijednostima, potrebno je to sve skupa prepisati na odgovarajući način u za to 


odredjene formulare. Na temelju.tih formulara buše se kartice i programiranje je završeno. 


Postoje tri grupe podataka: . j : 
a) Podaci okonfiguraciji odredjuju medjusobne veze izmedju pojedinih elemenata i spe- 
cificiraju željene operacije. i : Fo» 


b) Podaci o parametrima nadopunjuju podatke o konfiguraciji s numeričkim vrijednostima. 


c) Podaci za generatore funkcija definiraju ulazno-izlazne odnose za generatore funkcija. 


Podaci se mogu unijeti u računalo i direktno pomoću tastature na centralnoj jedinici. 
Taj način unošenja podataka bitno je sporiji nego pomoću bušenih kartica, pa se u praksi izbje- 
gava. j . E g 
Za vrijeme rada digitalnog računala, operator (odnosno onaj čiji se program upravo raču- 
na) mora se nalaziti kod centralne jedinice računala da bi mogao u odredjenom trenutku dati preko 
tastature računalu nove podatke neophodne za nastavak računanja. Osim toga, odvijanje programa 
može se zaustaviti na svakom mjestu da bi se npr. dodale nove instrukcije ili da bi se promijeni- 
li neki podaci koji su ranije uneseni bilo pomoću bušenih kartica, bilo putem tastature (sl. 1). 


Rezultat računanja bit će odštampan na Štampaču centralne jedinice u obliku numeričkih 
podataka,.a grafički prikaz odabrane veličine nacrtat će crtač krivulja (IBM 1627) na formatu 
10" x 10" (= 25,4x25,4 cm). . + : 


rogram_1130_CSMP ie specijalno za rješavanje općih diferencijalnih jednadž- 


P redvidjen 
bi pomoću posebno blokovsko ezika. Pro j sastavljen da se postavljeni zadatak rje€=- 
šava pomoću posebnog algoritma za sortiranje koji odredjuje redoslijed računanja numeričke in- 


šava pomoću posebnog algoritma za sortiranje koji ocreejnje rece 


tegracije. 


Buđuću da je 1130 CSMP organiziran kao "paralelni" jezik, nije od važnosti redoslijed 


unašanja podataka unutar jedne grupe podataka. (NE smiju se izmiješati npr. podaci o konfigura- 











ćiji s podacima o parametrima, jer su to dvije različite grupe.) Algoritam zazsortiranje ispitu- 
je da su sve vrijednosti u ređu, da li postoje sve potrebne ulazne varijable, da li je blokov- 


ska shema uredno sastavljena itd. Ovo ispitivanje traje nekoliko sekundi, a provodi se nakon sva- 
ke promjene konfiguracije simulacionog blokovskog dijagrama. Ukoliko nešto nije u redu, na štam- 
paču centralne jedinice ispisuje se odgovarajuća poruka. 


Položaj svakog elementa u blokovskoj shemi odredjen je zapravo brojem koji odgovara ele- 
mentu ispred njega. To nužno vodi na blokovsku shemu u obliku zatvorene petlje. U jednom zadatku 
može se pojaviti i više petlji, ali se u svakoj mora nalaziti barem jedan element s memorijom 





On 


preklopke on 

ZA elonnonanunornneamniie 

Ez [o]D]L21[51(41(51[61[2](51(21f0)(0](12)(12] [14 Ort 
DISK DISK ; PROGRAM| | IMM 
UNLOCKI| | READY START | | sTOP 

: 0 INT 

mol BGGODOODOOGBOOB e 

CHECK 2 Q .SToP | | RESET 
+1[-1f) 

i] legs] SGI CILICIJE Fogg] Frege 
SELECT | | CHECK IAR LOAD | 
=, 4 

G J 


po LII GICIOGDOCOOCI 


SI. 1 





(npr. integrator). Najčešća pogreška kod sastavljanja blokovske sheme je tzv. algebarska petlja. 
koja ne sadrži memorijski gI€Ement: OS LJ) € 7 


potrebno up I elemente V i Y koji sadrže mehanizam iteracije da bi se prekinula algebar- 
ska petlja. i 


Integracija je aproksimirana numeričkom metodom poznatom pod nazivom modoficirani Eule- 
rov postupak ili postupak drugog reda Runge-Kutta. Jednostavnost ovog postupka dolazi do izražaja 
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u slučaju sastavljanja specijalnih elemenata koji sadrže memoriju. U programu se primjenjuje cen- 
tralizirana integracija. Skup diferencijalnih jednadžbi obradjuje s€ Kao vektorska jednadžba. U 
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programu CSMP računa“se preko potprograma izlaz iz svih elemenata specificiranih u konftgu- 
raciji _i na taj se način dobiva _jeđan vektor. Ovaj se iznos vektora tada upotrebljava za izraču- 
žavalić Havt-vrijeđnosti VEREori varijabli stanja, što u stvari predstavlja izlaz iz sgih inte= 
gratora. 





Treba napomenuti da je jedinica vremena (sekunda, minuta, sat, dan ili dr.) implicitna u 
definiciji modela, a nije funkcija programa 1130 CSMP. Tokom simulacije ne smije se mijenjati je- 
dinica vremena. Sva vremena potrebna pri proračunu moraju biti izražena u istim jedinicama. 


Kao izvorni jezik upotrebljen je FORTRAN, a svi programi rade unutar 1130 MONITOR SISTEMA. 


Simulacija za koju je potrebno tniđesetak elemenata, uz 200 integracionih intervala traje 
2 do 4 minute što ovisi o vrsti upotrebljenih €lemenata. Promjena“npr:. jednog parametra zahtijeva 
manje od jedne minute, pa se tijekom jednog sata može (teoretski!) provesti 10 do 15 simulacija. 


PODRUČJE PRIMJENE 


Područje primjene obuhvaća inženjersko projektiranje i analizu, fizikalne probleme, bio- 
međicinska istraživanja i dr. 4 


Inženjer s područja regulacione tehnike može pomoću ovog programa i digitalnog računala 
provjeriti nekoliko različitih varijanti kod projektiranja nekog regulacijskog sustava. Može se 
provesti proračun servo-sistema_za npr. radarsku antenu za praćenje aviona_isto tako, kao i pro- 
račun postrojenja za valjanje aluminijskog lima ili nekog kemijskog reaktora. 


Fiziolog u medicinskom laboratoriju može npr. modelirati kardiovaskularni sustav. Simula- 
cijom mogu se provjeriti istraživanja o pretpostavljenom mehanizmu kojim je predstavljen promatra- 
ni fenomen. Odziv modela na simulirani vanjski podražaj omogućuje. usporedjivanje s kliničkim eks- 
perimentima. Efekti kirurških zahvata kao i bolesti često se puta mogu ispitati, provjeriti i na- 
dopuniti upravo pomoću modela. 


1130 CSMP ima široku primjenu u fizikalnim znanostima i nastavi. Studenti farmacije mogu 
na primjer, simulirati kemijski proces koji bi inaće bio neprikladan za neposredno proučavanje. U 
nastavi strojarstva, mogu se na taj način proučavati efekti prigušenja, zazora, zračnosti i slično 


= ll = 


promatranih mehaničkih naprava. Studenti elektrotehnike mogu provjeravati ponašanje sustava sa za- 
tvorenom petljom u koje spadaju i svi regulacijski sustavi. 


Sam program 1130 CSMP ima i pedagošku vrijednost, jer omogućuje stađenta da se više usre- 


dotoći na istraživani problem nego na mehanizam kojim se sustav simulira. 
činih Krako zadana: 
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OSNOVNA KONFIGURACIJA DIGITALNOG RAČUNALA 


Minimalna konfiguracija digitalnog računala IBM 1130 potrebna da se može upottijebiti 
program 1130 CSMP je slijedeća: 


1131 centralna procesna jedinica, Model 2B (memorija 8 k i disk) 








1442 čitač - bušač kartica 


1627 crtač krivulja. (Crtač krivulja nije neophodan, ali je vrlo koristan za ovu vrstu računanja). 


2. ELEMENTI ZA SASTAVLJANJE BLOKOVSKE SHEME 


U tablici 1 dani su elementi koji se mogu upotrijebiti prilikom sastavljanja blokovske 
sheme dinamičkog sustava kojeg se želi modelirati na digitalnom računalu IBM 1130. Postoji i pro- 
širena verzija programa CSMP za velike sisteme IBM 360 koja se razlikuje od 1130 CSMP u dodatnim 
algebarskim i logičkim mogućnostima digitalnog programiranja uz znatno veću fleksibilnost rada. 


Za linearne sustave često je puta dovoljno upotrijebiti samo prvih deset elemenata iz 
tablice 1. Ostali elementi služe uglavnom za simuliranje impulsnih i nelinearnih sustava, Za spe- 
cijalne elemente, kojih može biti najviše pet, mora se složiti posebni podprogram, 
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Generator funkcije F 





Drugi korijen H e e; e,=/e, 
Generator slučajnih J e. Generator slučajnih brojeva 
brojeva A izmedju + 1 





Apsolutna vrijednost M e 
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Ograničavač L —[1 >— e; 


Odstranjivanje nega- 
tivnih vrijednosti 


Odstranjivanje pozi- 
tivnih vrijednosti 


Relej 


Generator impulsa 


Jedinično kašnjenje 
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Generator impulsa kada je 2420 
e e; Perioda impulsa: Pi 
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Praznina 


Element za implicitne 
operacije 


Držanje nultog reda 


Specijalni elementi 








e, + £(e.) 
Upotrebljava se samo u spoju 
s elementom Y 


Logički element koji se upotreb- 
ljava za implicitne operacije 


= > 
e, e, za e, 0 


e, nepromijenjeno za e, š (4. 


Sadržaj odredjen odgovarajućim 
potprogramom 


imati svoj broj, tj. 
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Oznaka n označuje broj pojedinog elementa u cijeloj blokovskoj shemi. Svaki se element 
u blokovskoj shemi mora proizvoljnim redoslijedom označiti jednim brojem od 1 do 75 (jer je 
maksimalni broj elemenata koji se može upotrijebiti kod jednog računanja), ali svaki element mora 
dva ili više elemenata ne mogu imati jednake brojeve. 


kao nezavisne varijable rezerviran je broj 76. 


OZNAKA 
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ENGLESKI NAZIV 


BANG-BANG 

DEAD SPACE 
FUNCTION GENERATOR 
GAIN 

HALF POWER 
INTEGRATOR 

JITTER 

CONSTANT 

LIMITER 

MAGNITUDE 

NEGATIVE CLIPPER 
OFFSET 

POSITIVE CLIPPER 
QUIT 

RELAY 

TIME PULSE GENERATOR 
UNIT DELAY 


U tablici 2 dani su engleski nazivi za pojedine elemente. 


SLOBODNI PRIJEVOD NA HRVATSKI 


dvopoložajni element 

mrtva zona 

generator funkcije 

pojačalo 

drugt korijen 

integrator 

generator slučajnih brojeva 

konstanta 

ograničavač 

apsolutna vrijednost 

odstranjivanje negativnih vrijednosti 
pribrajanje konstante 

odstranjivanje pozitivnih vrijednosti 
gotovo 

relej 

generator impulsa 


jedinično kašnjenje 


Za definiranje vremena 
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V VACUOUS praznina 
." WEIGHTED SUMMER zbrajalo s množilom 
. x MULTIPLIER množilo 
Y WYE element za impliiitne operacije 
Z ZERO ORDER HOLD držanje nultog reda 
+ SUMMER zbrajalo 
/ DIVIDER djelilo 
. = SIGN INVERTER invertor 
l=5 SPECIAL specijalni elementi 


Kod mnogih elemenata važan je redoslijed ulaznih varijabli. Za elemente s jednim ulazom, 
ulaz se specificira kao "ulaz 1". Za elemente označene s 0, R, Y, Zi / ulaze treba uzeti kao što 
je navedeno u tablici 1. Kad bi se npr. zamijenili ulazi kod sklopa za dijeljenje dobili bi se 
recipročna vrijednost od tražene. Općenito se ulaz 3 ne smije upotrijebiti ako nisu upotrebljeni 
ulazi 1i2,a ulaz 2 se ne smije upotrijebiti ako nije upotrebljen ulaz 1. Izuzetak čine elemen- 
ti označeni s I, + iW. Kod njih se može ući na bilo koji ulaz. 


Kod integratora može se doduše ući na bilo koji ulaz (ako npr. postoji samo jedna ulazna 
veličina), medjutim treba zapaziti da prvi ulaz ulazi direktno u integrator, dok drugi i treći 
idu preko kruga koji označava množenje s parametrom 2, odnosno parametrom 3. To znači da je u slu- 
čaju, kađ je potrebno izlaz iz prethodnog elementa samo integrirati, najprikladnije ući na ulaz 1, 
jer bi za slučaj ulaska na ulaz 2 ili 3, trebalo još definirati iznos parametra P2, odnosno P3 
(koji bi u tom slučaju bili jednaki jedinici), i upisati ga u listu podataka o parametrima. i po- 
čČetnim uvjetima. 


Ako na ulaz integratora dolazi samo jedna veličina koju treba prethodno pomnožiti s ne- 
kom konstantom, onda treba ući na ulaz 2 ili na ulaz 3. Moglo bi se ući ina ulaz 1, ali bi prije 
njega trebalo upotrijebiti još jedan element tipa G (pojačalo) pomoću kojeg bi se provelo množe- 
nje. Medjutim u tom slučaju potrcšen je jedan element više nego što je potrebno, što treba izbje- 
gavati. , 


= 21 ša 


Treba naglasiti da parametri (prema tablici 1) pripadaju elementu i da ih ne treba ozna- 
čavati svojim brojem (jer su već odredjeni brojem elementa i brojem ulaza). Isto tako ne treba ih 
unositi u listu podataka o konfiguraciji, nego samo u listu podataka o parametrima i početnim uv- 
jetima. 


Kod elemenata W i + vidi se iz pripadne jednadžbe da je svejedno koji se ulaz upotrijebi 
kad postoji samo jedna (ili dvije) ulazne veličine. Neupotrebljeni ulazi i nespecificirani para- 
metri automatski se postavljaju na nulu. Razlika izmedju W i + je u tome, da u + ulazne veličine 
ulaze direktno, a kod elementa W moguće je pomoću parametara Pi, P2 i P3 pomnožiti svaki ulaz pro- 
izvoljnom konstantom. Element za zbrajanje (oznaka +) je jedini element kod kojeg je dozvoljen ne- 
gativni predznak kod specificiranja ulaza u listi podataka o konfiguraciji. 


Kod elemenata koji imaju memoriju (označeni s I, U, Vi Z), prvi je parametar upotrebljen 
za postavljanje početnih uvjeta, odnosno vrijednosti izlaza u t=0. I početni uvjeti koji nisu eks- 
plicitno specificirani, automatski se postavljaju na nulu. 


Pojedini element može imati najviše tri ulaza, ali svi imaju samo jedan izlaz, osim ele- 
menta Q koji nema izlaza. Jednom elementu mogu biti pridružena najviše tri parametra. Cijeli niz 
elemenata nema svoje parametre, a neki nemaju niti jedan ulaz (elementi SJ, Ki V). 


Jedini element kod kojeg je moguće bez posebnog sklopa (invertora) promijeniti predznak 
bilo koje ulazne veličine, je element za zbrajanje (oznaka +). Želi li se npr. dobiti razlika iz- 
medju izlaza iz nekog elementa označenog npr. s brojem 12 i izlaza iz elementa označenog s brojem 
26, onda se to na listu CONFIGURATION DATA (podaci o konfiguraciji) upisuje na slijedeći način: 


BLOK TIP ULAZ 1 ULAZ 2 : ULAZ 3 
16 + 12 =26 


Implicitne operacije 


Elementi označeni s Y i V upotrebljavaju se za prekidanje zatvorene petlje kad se prili- 
kom simuliranja zahtijeva neka implicitna operacija oblika _ y = f(x, y). Za ovakve odnose primi- 
jenjuje se kod računanja na digitalnom računalu postupak iteracije. Za y postavi se neka procije- 
njena početna vrijednost. Računanjem desne strane jednadžbe može se dobiti bolja procjena veličine 
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y. Ovaj se postupak ponavlja dok se ne dobije veličina y jednaka prethodno izračunatoj. Budući da 
se u općem slučaju i varijabla x mijenja tijekom računanja cijelog zadatka, postupak iteracije 
mora se ponoviti iza svakog intervala numeričke integracije. 


Procijenjena početna vrijednost veličine y predstavlja početni uvjet elementa V. Tije- 
kom daljnjeg računanja element V daje automatski daljnje procjene veličine y. Pridruženi element 
Y odredjuje da li je iterativnim postupkom razlika izmedju dviju procijenjenih vrijednosti veli- 
čine y postala manja od neke unaprijed zadane vrijednosti. Ako nije, onda nova vrijednost procje- 
ne mijenja izlaz iz elementa V i postupak iteracije se ponavlja. Nakon što je razlika rezultata 
dviju iteracija postala manja od željenog iznosa na izlazu iz elementa Y, dobiva se varijabla y 
i računalo nastavlja proračun u skladu s ostalim dijelom biokovske sheme. 


Kako uključiti elemente Y i V u blokovsku shemu na slici 2 prikazuje slika 3. Na slici 2 
prikazan je dio simulacione blokovske sheme u kojem treba obaviti računanje neke eksplicitne 
funkcije y = f(x, y). Kao što se vidi, pojavila se 
petlja jedinične povratne veze. Ovu petlju treba pre- 
kinuti na izlazu iz f(x, y) i prema slici 3 ubaciti 
elemente V i Y. Element V koji prema definiciji nema 
ulaza, služi samo kao spremnik probnih vrijednosti ve- 
ličine y. Izlaz iz V predstavlja ulaz u dio simulacio- 
nog dijagrama za računanje funkcije y = f(x, y). Ovaj 
dio simulacionog dijagrama može sadržavati i nekoliko 
različitih elemenata. Izlaz iz V ujedno se vodi na 
ulaz 2 elementa Y. Ulaz 1 u element Y je prekinuti iz- 
laz iz £(x, y). 





Prvim parametrom elementa Y odredjuje se dopušteni iznos razlike rezultata dviju uzas- 
topnih iteracija. Drugim parametrom odredjuje se iznos promjene nakon svake iteracije. To je u 
stvari neki faktor prigušenja. Element Y sadrži u sebi potprogram za izračunavanje iznosa 


1.0 - —>—-> 
Yn-1 


i usporedjivanje s iznosom prvog parametra. Kad orvi parametaxz postane veći od izračunate vrijed- 
nosti, zaustavija se daljnji postupak iteracije i računa se nova vrijednost veličine y prema re- 
laciji : 
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Parametar P2 može imati 
vrijednost izmedju nule i jedinice. 
Ako se implicitna operacija dobro 
ponaša i ako se konvergencija brzo 
postigne, može se postaviti P. jed- 
nako nuli Kad postupak iteracije 
teži oscilacijama preporuča se po- 
staviti P2 =0.5. 


Vrijednost parametra P2 
potrebno je odrediti eksperimentom. 
Ukoliko postupak iteracije traje 
suviše dugo, računanje se može pre- 
kinuti pomoću preklopke 0, da bi se 
u računalo unijela nova veća vri- 
jednost parametra P2. 


Isto tako, da bi se povećala efikasnost procesa iteracije, potrebno je da element V (od- 
nosno elementi V, ako ih ima više) budu označeni s najvišim brojevima koji označuju broj elementa. 


Impulsni sustavi 


Za simuliranje impulsnih sustava koji za proračun upotrebljavaju jednadžbe diferencija i 


Z-transformaciju, razvijeni su elementi označeni s T, Zi. 


Generator impulsa (element T) proizvodi impulse čija je perioda odredjena parametrom Pi: 
Ovaj će element davati na svom izlazu slijed impulsa tako dugo dok mu je ulaz veći ili jednak nu- 
li. Ovim generatorom impulsa najlakše se upravlja pomoću bloka 76. Blok 76 predstavlja prema de- 
finiciji nezavisnu varijablu, vrijeme. Izlaz iz generatora impulsa obično se upotrebljava kao 
ulaz 2 (triger) u element Z. Za vrijeme pozitivnog impulsa iz elementa T, element Z uzima uzorak 
svoje ulazne vrijednosti i zadržava ga na svojem izlazu do slijedećeg impulsa iz T. Ova kombina- 
cija elemenata T i Z predstavlja impulsni element u blokovskoj shemi impulsnog sustava. Jedan 


generator može okidati i više elemenata Z. 


pz 24 - 


Ako se želi simulirati prijenosna funkcija višeg reda, potrebno je izmedju pojedinih Z 
elemenata postaviti elemente U (jedinično kašnjenje) da bi se dobili izrazi Z-i, Z-2 itd. Element 
U zadržava ulaznu vrijednost za element Z koji je spojen na njegov izlaz. Kad ne bi postojao ele- 
ment U izmedju dva Z elementa, izgubila bi se izlazna vrijednost prethodnog Z elementa, jer raču- 
nalo računa izlaze iz pojedinih elemenata po redoslijedu toka informacije. Element U ima kašnje- 
nje jednako polovici integracionog intervala. : 


Izmedju vremenskog kašnjenja od KZ elementa U i intervala uzimanja uzoraka odredjenog 


periodom impulsa iz elementa T nema medjusobnog odnosa. 


Generatori funkcija 


Generator funkcije (element F) služi za modeliranje proizvoljnih funkcija oblika y = f(x) 
koje se mogu aproksimirati s deset pravocrtnih segmenata. Maksimalna i minimalna vrijednost ulazne 
varijable x u promatranom intervalu, mora se zadati u listi podataka o parametrima i početnim 
uvjetima, a vrijednost funkcije na kraju svakog od 10 odsječaka mora se definirati na listi poda- 
taka o generatorima funkcija (lista: FUNCTION GENERATOR DATA). Maksimalno se može kod jednog raču- 
nanja upotrijebiti tri generatora funkcija. 


Specijalni elementi mogu se definirati prema želji i potrebama. To može biti element koji 
će npr. predstavljati histerezu, funkciju sinusa, kosinusa, eksponencijalnu funkciju itd. Za ovak- 
ve seecijalne elemente ili već postoje posebni programi, ili ih treba tek napisati u FORTRANU. 
Maksimalni broj specijalnih elemenata u jednom računanju je pet. 


3. PRIPREMA PODATAKA 


3.1. CRTANJE BLOKOVSKE SHEME 


Postavljeni zadatak treba najprije nacrtati u obliku blokovske sheme upotrebljavajući 
elemente prema tablici 1. Pritom treba paziti na ograničenja i pravila spomenuta u prethodnom po- 
glavlju. 


Prikladno je na samu blokovsku shemu napisati osim rednog broja elementa. i njegovu ozna- 
ku, a pored kruga koji označava pojedini parametar i njegovu numeričku vrijednost. Zbog lakšeg 
snalaženja u blokovskoj shemi, može se iza svakog elementa staviti oznaka ili naziv veličine koja 
se vodi u slijedeći element (vidi primjere). 


3.2. PODACI O KONFIGURACIJI 


Na listu CONFIGURATION DATA (podacet o konfiguraeiji)* sl.4, upisuju se podaci o elementi- 
ma blokovske sheme. Svaki element upisuje se u svoj redak, a definiran je slijedećim podacima: 


a) naziv izlaza Cc) oznaka 
b) broj elementa d) broj elementa koji se nalazi ispred 


Kosim slovima pisani su slobodni prijevodi engleskih izraza . 
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a) Naziv izlaza (OUTPUT NAME) može se i ne mora napisati. Taj se naziv obično upisuje 
samo za one izlaze koji su interesantniji za sam proračun, a za ostale se ostavlja prazan prostor. 
Za ovaj naziv rezervirano je prvih 16 znakovnih mjesta na formularu, što ujedno znači da se karti- 
ca mora bušiti na tim poljima. 


Sedamnaesto i osamnaesto mjesto mora ostati prazno. 


b) Broj elementa (BLOCK) treba upisati u polja 19 i 20. Brojeve elemenata od 1 do 9, tj. 
one koji imaju samo jednu znamenku, treba upisati u polje 20. Za broj elementa mogu se upotrijebiti 
samo brojevi od 1 do 75. Broj 76 rezerviran je za element kojim se može definirati vrijeme kao ne- 
zavisna varijabla. 


Polja 21 - 29 ostaju prazna. 


c) Oznaku elementa (TYPE) treba upisati u polje 30. 
Polja 31 - 37 ostaju prazna. 


d) Ulazi u element (INPUT 1, INPUT 2, INPUT 3) odredjeni su brojem elementa koji se na- 
lazi ispred. Ovi se brojevi upisuju u polja 38 - 40, 48 - 50 i 58 - 60. Brojeve treba upisati tako 
da budu poravnani s desne strane. To znači da će se brojevi elemenata od 1 do 9 upisivati u polja 
40, 50 ili 60, a dvoznamenkasti brojevi u polja 39 - 40, 49 - 50 odnosno 59 - 60. Treće polje (38, 
48 i 58) rezervirano je za eventualni znak minusa koji se može pojaviti samo kod elementa za zbra- 
janje (oznaka +). 


Polja 61 - 72 ostaju prazna. 


Bušene kartice s podacima o konfiguraciji mogu biti složene u bilo kakvom redoslijedu, 
ali se moraju nalaziti zajedno i ne smiju se pomiješati s drugim bušenim karticama. 
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3.3. PODACI O _ PARAMETRIMA 


Na listu INITIAL CONDITIONS AND PARAMETER DATA (početni uvjeti i podaci o parametrima) 
sl. 5, upisuju se numeričke vrijednosti početnih uvjeta i parametara koji su pridruženi pojedinom 
elementu. U ovu listu upisuju se samo oni elementi kod kojih su definirani početni uvjeti i para- 
metri. Za takve elemente upisuju se slijedeći podaci: 


a) naziv početnog uvjeta ili parametra 
b) broj elementa 

C) početni uvjet ili prvi parametar 

d) drugi parametar 

e) treći parametar 


a) Naziv početnog uvjeta ili parametra (IC/PAR NAME) upisuje se u prvih 16 polja. U ova 
se polja može upisati i naziv elementa kojem pripadaju podaci. Polja se mogu ostaviti i prazna, 
jer ovi podaci ne sudjeluju u računanju, nego služe samo za lakše snalaženje. 


Polja 17 - 18 ostaju prazna. 


b) Broj elementa (BLOCK) upisuje se u polja 19 i 20. Za način upisivanja broja u polja 
vrijede ista pravila kao i za upisivanje u listu podataka o konfiguraciji. 


Polja 21 - 25 ostaju prazna. 


€) U polja 26 - 35 treba upisati numeričku vrijednost početnog uvjeta ili prvog pakame- 
tra (I/C PAR 1). Početni uvjet mogu imati samo elementi s memorijom (označeni s I, U, Vi Z). Za 
elemente koji nemaju početne uvjete u ova se polja upisuje iznos prvog parametra. Zapravo do za- 
bune ne može niti doći, jer su i početni uvjeti označeni s Pi kao i prvi parametar. 


Polja 36 - 40 ostaju prazna. 

d) Podatke za drugi parametar (PAR 2) treba upisati u polja 41 - 50, a za treći parame- 
tar (PAR 3) u polja 56 - 65. 

Polja 51 - 55 i 66 - 72 ostaju prazna. 
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I ovdje vrijedi pravilo da podatke za svaki element treba pisati u svoj redak. Numeričke 
vrijednosti parametara treba upisati kao realne konstante. Nužno je upisati decimalnu točku. Broj 
koji označuje vrijednost parametra može se upisati bilo gdje unutar predvidjenih polja (26 - 35, 
41 - 50 odnosno 56 - 65), tj. nije ga potrebno poravnavati npr. uz desni rub. Neka npr. prvi pa- 
rametar ima iznos 2, onda treba u bilo koja dva polja izmedju dvadesetšestog i tridesetpetog upi- 
sati 2. , a računalo će samo smjestiti taj broj u polje 30, decimalnu točku u polje 31, a u polja 
31 - 35 upisat će četiri nule. Ako se podaci unose u računalo pomoću bušenih kartica, onda je do- 
pošten i E format. 


Bušene kartice s podacima o parametrima mogu biti poredane u bilo kojem redoslijedu, ali 
se moraju nalaziti zajedno. 


3.4. PODACI O GENERATORIMA FUNKCIJA 


Na listu FUNCTION GENERATOR DATA (podaet za generatore funkeija) sl. 6, upisuju se iznosi 
funkcije na krajevima pojedinih intervala i broj elementa. Za svaki generator funkcije (kojih može 
biti najviše tri) potrebne su tri bušene kartice, jer se podaci za svaki generator funkcije upisu- 
ju u tri reda. 


U prvi redak upisuje se: u polja 19 i 20 broj elementa (BLOCK) na isti način kao i na li- 
stu podataka o konfiguraciji; u polja 26 - 35 početna vrijednost funkcije, što je ujedno, jednako 
vrijednosti funkcije za x = P2 (tj. treba upisati iznos parametra P2); u polja 41 - 50 iznos 
funkcije na kraju prvog, u polja 56 - 65 iznos funkcije na kraju drugog i u polja 71 - 80 iznos 
funkcije na kraju trećeg odsječka. 


U drugi redak (a prema tome: na drugu bušenu karticu) upisuju se u polja 26 - 35, 41 - 50 
56 - 65 i 71 - 80 iznosi funkcije na kraju četvrtog, petog, šestog i sedmog odsječka. 


U treći redak i to u polja 26 - 35, 41 - 50 i 56 - 65 upisuju se iznosi funkcije iza po- 
sljednja tri odsječka. Iznos funkcije iza zadnjeg odsječka ujedno je i iznos funkcije za x =P1, 
gdje je Pi, iznos prvog parametra. : 


Bušene kartice s ovim podacima moraju biti složene po redu i ne smiju se izmiješati ni 
medjusobno ni s ostalim bušenim karticama. 
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Ukoliko se ne raspolaže s posebnim formularima za upisivanje podataka (sl. 3,415), po- 
daci se mogu upisati i na obične formulare za programiranje u FORTRANU (sl. 7), ali se mora paziti 
da se podaci upišu točno u odredjena polja. 


3.5. PRIPREMA BUŠENIH KARTICA 


Prije nego što se pristupi radu na računalu potrebno je pripremiti još nekoliko bušenih 
kartica. Kompletni paket bušenih kartica koje treba dati računalu sadrži: 


(a) startna kartica (COLD START CARD) 

b) XEQ kartica 

c) dvije LOCAL kartice- 

d) kartice s podacima o konfiguraciji 

e) jedna prazna kartica š 

f) kartice s podacima o početnim uvjetima i parametrima 
g) jedna prazna kartica 

h) kartice s podacima za generatore funkcija 

i) nekoliko praznih kartica 


a) Startna kartica isporučuje se zajedno s diskom za program 1130 CSMP. Na slici 8 vidi 
se kako je izbušena ova kartica. 


b). XEQ kartica (sl. 9) daje nalog da se program 1130 CSMP prebaci iz diska u magnetsku 
memoriju i da se na nju prebaci upravljanje. Na ovu karticu treba izbušiti slijedeće: na prva dva 
mjesta dolazi kosa crta (/); treće mjesto je prazno; zatim dolazi XEQ; sedmo mjesto je prazno, a 
iza njega dolazi CSMP; polja od 12 - 16 su prazna, a na polje 17 treba upisati broj 2; polja od 
18 - 72 ostaju prazna. : 


€) LOCAL je kratica za LOADED ON CALL. Pomoću ovih kartica pozivaju se s diska i upisuju 
u magnetsku memoriju potprogrami koji se najčešće upotrebljavaju tijekom računanja. Svi ostali 
Potprogrami pozivaju se prema potrebi s diska. 


Re 34 - 
Na prvu LOCAL karticu (sl. 10) treba u polja 1 - 69 izbušiti:  *LOCALCSMP,CSMO,CSMI, 
CSM2,CSM3,CSM4,CSM5,CSM6,CSM7,CSMBA,CSMII1, CSMI2., 
Polja 70 - 72 ostaju prazna. 
Ako disk za program 1130 CSMP ima program sastavljen za rad s crtačem krivulja 1627, na 
drugu LOCAL karticu (sl. 11) treba u polja 1 - 23 izbušiti:  *LOCALCSM13,CSMJA,CSM9B 


Polja 24 - 72 ostaju prazna. 
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Ukoliko digitalno računalo nema crtača krivulja 1627, na ovu drugu LOCAL karticu treba u 
polja 1 - 11 izbušiti:  xLOCALCSMI3 


Polja 12 - 72 ostaju prazna. 


li 
00000000000000000000000000000000000000000000000000000000000000000000000 
E ITET TIE TEKELEEEREEDEEELELEEDEUDE KH 
E EE E EEE EEE EEE EEEEEEEEEEEEEARREERKREEAREEE 
PETETEVI VETETETI VETETE TETE EEE ET TE ETE TETI ETE ETE TETE KEKE KEZELE EE EEK EKE EEE KEKE 
3333333333333333333333333333333333333339333333333333333333333333333333333339333 
mrve vrererrerrrrereverrerrere tet Tea rta tree z zeza ZEKE EEE EEE EEE EEAS 


$555m5555555555555555555555555555555555595555555555553555595595555533559339533%5 


 D 
u 
ra 
š 
x 
z 
IT) 
u 
> 


66666666666666666666666666666666666666666666666666666666666666666666666666GEEGGH 
EE VREREI UREREE REDE RERRA EE REKA RRRRARARAKERRAKERARERARERARERREKEREKERSEESEKEKEE 


88888g08886888888888868888886888888888B6B8BBBBBBBBEBB6BBBBBBBBBEBBBBHBEBBEBKAIKI 


3999999 99999 
12345617 6 70 16 1928 





LOCAL karticama označuju se podprogrami koji će ostati na disku i tek kad se ukaže potre- 
ba za njima onda ih se poziva. Ovo se radi zato da. se centralna memorija od 8 k ne bi nepotrebno 
nakrcala trenutno nepotrebnim podacima. Medjutim kako digitalno računalo IBM 1130 na ETF-u ima me- 
moriju od 32 k, štednja ove vrste postaje suvišna, pa prema tome nisu niti potrebne LOCAL kartice. 
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U tom slučaju i kartica XEQ nema u koloni 17 izbušenu brojku 2 (sl. 12), koja inače označava da 
iza kartice XEQ dolaze dvije LOCAL kartice. 
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l 
d) Kartice s podacima o konfiguraciji mogu biti poredane bilo kojim redoslijedom, ali se 
moraju nalaziti zajedno. 


f) Kartice s podacima o parametrima mogu biti takodjer poredane u proizvoljnom redosli- 
jedu. Važno je samo, da se sve nalaze zajedno. 
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h) Kartice s podacima za generatore funkcija moraju biti složene tako da se za svaki ge- 
nerator funkcije po tri kartice moraju nalaziti složene po redoslijedu kako su upisivani podaci u 
listu podataka za generatore funkcija. U protivnom pomiješali bi se iznosi funkcije iza pojedinih 
odsječaka. 
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e), g) Prazna kartica iza svake grupe podataka označuje kraj va Teoma ja te grupe podata- 
ka i ujedno odjeljuje jednu grupu od druge. 

Ukoliko se u programu za rješavanje postavljenog zadatka ne nalaze generatori funkcija, 
onda se stavke g) i h) mogu izostaviti. 
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i) Za jedno računanje dovoljna bi bila samo jedna prazna kartica .na kraju kompletnog pa- 
keta kartica. Medjutim rezultat prvog računanja redovno ne zadovoljava, pa treba :u računalo unije- 
ti promjenu nekog podatka da bi se dobio bolji rezultat. Budući da računalo uvijek pokuša dobiti 
podatak iz čitača kartica prije nego dopusti ulazak podataka pomoću tastature, mora se uvjeriti da 
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iza prazne kartice koja je označila kraj zadnje grupe kartica, nema više podataka na bušenim kar- 
ticama, Tek nakon što očita praznu karticu, moguće je unijeti novi podatak pomoću tastature. Zah- 
tijeva li proračun mnogo promjena podataka, može se dogoditi da je potrebno u čitač kartica dodati 
još praznih kartica, 
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Ukoliko se unaprijed znađe koje promjene treba učiniti nakon prvog računanja, mogu se no- 
vi podaci pripremiti na bušenim karticama. Novi podaci mogu se pomoću bušenih kartica pripremiti 
samo u vezi promjene podataka o konfiguraciji, početnih uvjeta i parametara te podataka o genera- 
torima funkcija. 


Ako se npr. nakon prvog računanja želi provesti još nekoliko računanja, ali s drugim una- 
prijed odabranim parametrima, potrebno je pripremiti snop kartica i to tako, da se iza svake kar- 
tice (ili grupe kartica) s novim podacima, iza kojih se traži ponovno računanje, nalazi po jedna 
prazna kartica. Očitavanje novih podataka o parametrima obavit će se nakon prebacivanja preklopke 
broj 2 u gornji položaj i pritiskivanja na tipku PROGRAM START. (Za promjenu podataka o konfigu- 
raciji treba prebaciti preklopku 1, a za promjenu podataka za generator funkcije, preklopku 3 u 
gornji položaj.) 


Detaljnija objašnjenja o grupi kartica poda), b) ic), kao i o njihovoj funkciji mogu 
se naći u [4]. 


4. RAD NA RAČUNALU 


Nakon što je kompletiran snop bušenih kartica može se pristupiti radu na računalu. Taj se 


rad .može podijeliti na nekoliko faza: 
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1) 
2) 
3) 
4) 
5) 
6) 
7) 
8) 


pripremna faza 

početna faza 

faza primanja podataka o konfiguraciji 

faza primanja podataka o parametrima 

faza primanja podataka o generatorima funkcija 
faza davanja uputa računalu 

faza računanja 

faza zahtjeva po želji 


PRIPREMNA FAZA 


CSMP, 
tipke RUN, IMM STOP i RESET. 


m. 


Kao prvo potrebno je u centralnu jedinicu računala staviti disk 2315 s programom 1130 
(Na ETF-u taj disk ima broj 5002). Zatim treba pritisnuti tipku FILE SWITCH i na konzoli 
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Prebacivanje programa iz diska u magnetsku memoriju i priprema za- početak prijema podata- 
ka traje cca 3,5 minute. Za to vrijeme linijski štampač ispisuje (uz pretpostavku da su na početak 
snopa bušenih kartica stavljene prve četiri bušene kartice prema opisu u 3.5./: 


PAGE 1 
// 0B 5002 


LOG DRIVE CART SPEC CART AVAIL  PHY DRIVE 
0000 5002 5002 0000 


V2 MO7 ACTUAL 32K  CONFIG 32K 
// XEQ CSMP 2 
*LOCALCSMP 9 CSMO »CSMI »CSM2 »CSM3 #CSM4 9 CSM5 9 CSM6 #CSM7 »CSMBA s CSMI 1 #CSMI2» 


#LOCALCSM13 »CSM9A » CSM9B 





4.2. POČETNA FAZA 


Odmah nakon što je računalo obavilo sve potrebne predradnje, program počinje sa štampa- 
njem naslova 


CONTINUOUS SYSTEM MODELING PROGRAM 
A DIGITAL ANALOG SIMULATOR PROGRAM FOR THE IBM 1130 


(Program za modeliranje kontinuiranih sustava 
Program za analogno-digitalnu simulaciju za digitalno računalo IBM 1130) 


Treba naglasiti da se ovaj i sav ostali tekst kao i podaci štampaju na štampaču koji se 
nalazi na centralnoj jedinici. Linijski Štampač mora takodjer biti uključen, ali se na njemu ne 
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štampa više ništa. (Štampao se je tekst samo za vrijeme očitavanja prvih četiruju, odnosno dviju, 
bušenih kartica.) j 


Nakon naslova odštampaju se:četiri upute: 


INSTRUCTI ONAL COMMENTS MAY BE SUPPRESSED AT ANY TIME BY TURNING ON SWITCH 10 
(Štampanje uputa može se prekinuti u svakom trenutku uključivanjem preklopke 10) 


Postavljanjem preklopke 10 u gornji položaj može se izbjeći štampanje nekih uputa koje 
nisu neophodne i koje uglavnom služe kao podsjetnik koju preklopku treba uključiti, da bi se pro- 
gram odvijao na željeni način. Ove se upute može lako naučiti napamet, pa nije potrebno da ih ra- 
čunalo svaki put ispisuje, jer se time nepotrebno gubi vrijeme. I ova se uputa neće odštampati, 
ako se preklopka 10 uključi na vrijeme, tj. prije nego li je počelo njezino štampanje. Uključi li 
se preklopka 10 za vrijeme štampanja poruke, ona će se odštampati do kraja, ali se slijedeće upu- 
te ne će štampati za cijelo vrijeme dok je preklopka 10 uključena. Preklopka 10 može se uključiti 
u bilo kojem trenutku. 


Ukoliko se preklopku 10 ne uključi, iza prve upute, bez zadržavanja štampa se druga uputa: 


TURN ON SWITCH 1 TO ENTER OR MODIFY CONFIGURATION STATEMENTS VIA THE KEYBOARD 

(Uključi preklopku 1 ako treba podatke o konfiguraciji unijeti ili promijeniti pomođu tastature) 
Ako se podaci o konfiguraciji unose u računalo samo pomoću bušenih kartica i čitača kar- 

tica, preklopka 1 mora biti isključena. Unose li se ti podaci, ili samo dio tih podataka, pomoću 


tastature, treba uključiti preklopku 1. U tom slučaju računalo će se zaustaviti prije štampanja 
naslova za slijedeću grupu podataka i omogućiti unošenje podataka preko tastature. 


I ova se poruka ne će odštampati ako je bila preklopka 10 uključena; ako nije bila uk- 
ljučena, odštampat će se bez zadržavanja treća uputa: 


TURN ON SWITCH 2 TO ENTER OR MODIFY INITIAL CONDITIONS OR ELEMENT PARAMETERS VIA THE KEYBOARD 


(Uključi preklopku 2 ako treba početne uvjete ili podatke o parametrima unijeti ili promijeniti 
pomoću tastature) 


Postupak je isti kao i kod prethodne upute, tj. preklopka 2 mora biti isključena ako se 
svi podaci ove grupe unose preko bušenih kartica i čitača kartica, a uključena ako treba ove 


= 
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podatke ili dio podataka unijeti u računalo ili promijeniti ih pomoću tastature. Nije li uključena 
preklopka 10, odštampat će se i četvrta uputa: 


TURN ON SWITCH 3 TO ENTER OR MODIFY FUNCTION GENERATOR INTERCEPTS VIA THE KEYBOARD 


(Uključi preklopku 3 ako treba podatke o generatoru funkcije unijeti iLi promijeniti pomoću tasta- 
ture) 


Za položaj preklopke 3 vrijede ista pravila kao i za prethodna dva slučaja. I preklopka 3 
može se isto kao i preklopke 1 i 2 uključiti u bilo kojem trenutku. 


U slučaju kad su svi podaci uredno pripremljeni na bušenim karticama, za vrijeme ove faze 
nije potrebno izvršiti nikakve radnje. Naslov i sve četiri poruke odštampat će se najvećom mogućom 
brzinom i program će se nastaviti trećom fazom. 


4.3. FAZA PRIMANJA PODATAKA O KONFIGURACIJI - 


Ova faza počinje štampanjem naslova 


CONFIGURATION £PECIFICATION 


(Podaci o kcsxfiguraciji) 


Nakon toga ispisuju se nazivi pojedinih ko.ona +2ste CONFIGURATION DATA (podaci o konfiguraciji): 


OUTPUT NAME BLOCK TYPE INPUT 1 INPUT 2 INPUT 3 
(Naziv izlaza broj elementa oznaka ulaz 1 nlaas 2 ulaa €; 


Tek sada počinje očitavanje kartica s podacima o konfiguraciji koji se istovremeno štam- 
paju u odgovarajuće kolone. Da bi se uštedjelo na vrekcic4, zožše se isključiti štampanje ovih poda- 
taka uključivanjem preklopke 11. U tom slučaju očitavanje kartica ide znatno brže. Preklopka 11 
može se uključiti u bilo kojem trenutku. Ako se uključi za vrijeme štampanja ovih podataka, onda 








se slijedeći redak više neće odštampati. 


Računalo će očitavati kartice i spremati podatke u memoriju dok ne naidje na prvu praznu 
karticu koja označava kraj ove grupe podataka. Ukoliko je preklopka I bila uključena, računalo se 
zaustavlja spremno da preko tastature primi još neke podatke u vezi konfiguracije. 


Da bi se olakšalo unošenje podataka pomoću tastature, računalo će odštampati pet pari za- 
grada u koje treba prema ranije opisanim pravilima, upisati podatke o novom elementu ili promjenu 
podataka za već specificirani element. Zagrađe će se odštampati ovako: 


OUTPUT NAME BLOCK TYPE INPUT 1 INPUT 2 INPUT 3 
( ko) ( ) ( ) ( ) 


Nakon ovoga računalo stane i čeka dok se ne upišu potrebni podaci. Ako se slučajno pogri- 
ješi kod upisivanja podataka, bilo da se upiše krivi broj ili da se ga ne napiše unutar zagrada 
(što zahtijeva stanovitu vježbu u "gadjanju"), potrebno je ispravni podatak.napisati ponovno u no- 
vi redak. Na taj će način stari podatak biti automatski izbrisan. To će biti omogućeno pritiskom 
na tipku ERASE FIELD, jer će se tada glava za pisanje vratiti u početni položaj, a valjak će po- 
maknuti papir na novi redak. Ako se želi samo izbrisati npr. neki suvišno specificirani element, 
dovoljno je u novi redak napisati samo broj elementa, a ostale zagrade ostaviti prazne. 


Da bi se označio kraj unošenja podataka o konfiguraciji preko tastature, potrebno je is- 
ključiti preklopku 1 (vratiti je u donji položaj) prije nego li se u zagrade upišu podaci za zad- 
nji element. Tako dugo dok je preklopka 1 uključena (tj. dok se nalazi u gornjem položaju), raču- 
nalo očekuje daljnje upisivanje podataka preko tastature i automatski ispisuje zagrade za upisiva- 
nje novih podataka. 


Kad su uneseni svi podaci u vezi konfiguracije, potrebno je pritisnuti tipku EOF da bi se 
upisani podaci unijeli u program. (Ovdje se misli samo na podatke koji su uneseni preko tastature). 
Računalo tada analizira ispravnost unesenih podataka i sredjuje ih prema odredjenom algoritmu. 
Ukoliko ne pronadje nikakvu pogrešku, računalo prelazi na slijedeću fazu. Ako nešto nije bilo u 
redu upisano, računalo će ispisati vrstu pogreške i tražiti da se ona ispravi. 


Ispiše li računalo: 


THE SPECIFICATION FOR BLOCK NN IS INCORRECT 


(Specifikacija za element NN nije u redu) 
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to znači da za taj element nije dobro upisan njegov vlastiti broj ili broj nekog ulaza. Ova će se 
poruka ispisati ako se upisani broj ne nalazi izmedju 1 i 75, ili ako nije upisan uz desni kraj 
raspoloživih poija. Nakon ove poruke računalo ispiše pet pari zagrada unutar kojih treba sve po- 
datke za dotični element pravilno upisati pomoću tastature. Tek nakon što su podaci ispravno une- 
seni, računalo nastavlja rad. 


Ako se brojevi elemenata od 1 do 7, kod upisivanja ne poravnaju uz desnu zagradu, raču- 
nalo će ih shvatiti kao brojeve 10, 20, 30, ..., 70. Tek za npr. broj 8 ili 9 upisan uz lijevu za- 
gradu, ispisat će se gornja poruka, jer će to računalo shvatiti kao broj 80, odnosno 90, koji su 
veći od 75. 


ILLEGAL SYMBOL SPECIFIED FOR BLOCK NN 
(Za element NN upisana nepropisna oznaka) 
Ova će se poruka ispisati u slučaju kad se u polje 30 upiše oznaka elementa koja nije na- 


vedena u tablici 1. Nakon ove poruke računalo ispisuje pet pari zagrada u koje treba upisati is- 
pravne podatke. 


Da bi se ispravila npr. samo jedna pogreška u jednom retku, moraju se u novi redak ponov- 
no ispisati svi podaci tog retka, a ne samo podatak koji je pogrešno upisan. i 


PREVIOUS SPECIFICATION DELETED 


(Prethodna specifikacija je izbrisana) 


Ova će se poruka ispisati uvijek kad se podaci konfiguracije prethodnog elementa izbrišu, 
promijene ili još jednom napišu. 
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BLOCK NN, ONE OF THE INPUTS TO BLOCK MM, HAS NOT BEEN SPECIFIED 
(Element NN, jedan od ulaza u element MM nije naveden) 


Ova poruka znači da jedan od ulaza u element MM (tj. broj prethodnog elementa NN) nije 
definiran u podacima o konfiguraciji. Pogreška se sastoji u tome da je za ulaznu veličinu u ele- 
ment MM napisan broj elementa koji ne postoji u blokovskoj shemi (odnosno na listi podataka o 
konfiguraciji). Tada je potrebno ili ispraviti podatak za ulaznu veličinu u element MM, ili ispra- 
viti blokovsku shemu i tada ujedno u listu podataka o konfiguraciji dodati element koji nedostaje. 
Ova pogreška uglavnom nastaje kad se zabunom upiše (odnosno izbuši) krivi broj koji označava ulaz- 
nu veličinu, ili kad se zaboravi iz blok sheme prepisati u listu podataka o parametrima neki ele- 
ment (odnosno ako se koja kartica izgubi). 


Ova će se poruka odštampati umjesto podataka za element NN, a ispod nje odštampat će se 
pet pari zagrada u koje treba upisati ispravne podatke. 


1130 CSMP REQUIRES AT LAST 1 INTEGRATOR 
(Program 1130 CSMP zahtijeva barem 1 integrator) 


Budući da se u podacima o konfiguraciji mora nalaziti barem jedan integrator, pogreška 
se može ispraviti uključivanjem barem jednog integratora u blokovsku shemu. 


THE HAXIIUH GE. 25 UNIT BELAYS HAS LEEH EXCEEĐED 
(Upotrebljeno je više od 25 jediničnih kašnjenja) 


Moguće je upotrijebiti najviše 25 elemenata s oznakom U. Ako ih se upotrijebi više, po- 
trebno je preraditi blokovsku shemu da se njihov broj smanji. 
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THE MAXIMUM OF 3 FUNCTION GENERATORS HAS BEEM EXCEEDED 
(Upotrebljeno je više od tri generatora funkcija) 

Ova će se poruka odštampati ako je upotrebljeno više od tri generatora funkcija. Pogreška 
se otklanja smanjenjem broja generatora funkcija. Ako je neophodno raspolagati s više od tri gene- 


ratora funkcija, onda se mogu slijedeći generatori funkcija isprogramirati kao specijalni elementi 
(kojih može biti najviše pet). 


THE MAXIMUM OF 25 INTEGRATOSS HAS BEEN EXCEEDED 
(Upotrebljeno je više od 25 integratora) 


Budući da se ne može upotrijebiti više od 25 integratora, mora se preraditi blokovska 
shema da bi se poareška otklonila. 


SORT FAILURE - BLOCK NN 
(Pogreška razvrstavanja - element NN) 


Ako računalo nije sa svojim algoritmom za razvrstavanje u stanju odrediti redoslijed ra- 
čunanja, tj. ako sa ulaze u element NN dobiva druge podatke nego što su dani u podacima o konfigu- 
raciji, potrebno je promijeniti podatke o konfiguraciji da bi računalo moglo nastaviti s računa- 
njem. 


Treba napomenuti da i za vrijeme ove faze nije potrebno izvršiti nikakve radnje, ako su 
podaci na bušenim karticama uredno pripremljeni, pa se program bez zaustavljanja nastavlja na čet- 
vrtu fazu. 


4.4. FAZA PRIMANJA PODATAKA O PARAMETRIMA 


Ova faza počinje štampanjem naslova: 


INITIAL CONDITIONS AND PARAMETERS 
(Početni uvjeti i parametri) 


Nakon toga ispisuju se naslovi pojedinih kolona liste INITIAL CONDITIONS AND PARAMETER 
DATA (Početni uvjeti + podaet o parametrima): 


IC/PAR NAME BLOCK IC/PARI PAR2 PAR3 


Sada počinje očitavanje kartica s podacima o početnim uvjetima i parametrima koji se is- 
tovremeno štampaju u odgovarajuće kolone. Očitavanje kartica i spremanje podataka u magnetsku me- 
moriju nastavlja se sve dok čitač kartica ne naidje na prvu praznu karticu koja ujedno označava 
kraj ove grupe podataka. Računalo ne će ispisivati ove podatke ako se uključi preklopka 11. 


Nakon što su očitane sve kartice ove grupe podataka, računalo ispituje da li je uključe- 
na preklopka 2. Ukoliko nije, program se nastavlja ispitivanjem primljenih podataka. Pavovremenim 
uključivanjem preklopke 2 (tj. najkasnije za vrijeme štampanja posljednjeg retka s podacima o po- 
četnim uvjetima i parametrima), računalo se zaustavlja spremno da preko tastature primi još neke 
podatke u vezi početnih uvjeta i parametara. x 


Da bi se olakšalo unošenje i ovih podataka u odredjena polja, računalo će odštampati če- 
tiri para zagrada (u koje treba upisati nove podatke) ispod naslova pojedinih kolona, na slijede- 
ći način: : 


IC/PAR NAME BLOCK IC/PARI PAR2 PAR3 
) ) ( 


Nakon što se unesu podaci za željeni element, računalo ponovno odštampava zagrade u novi 
redak, i tako se nastavlja upisivanje podataka pomoću tastature. Pogriješi li se kod upisivanja 
podataka, treba pritisnuti tipku ERASE FIELD i ponovno ispravno upisati podatke za cijeli redak. 
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Da bi se označio kraj unošenja ovih podataka preko tastature, potrebno je isključiti pre- 
klopku 2 prije nego se upišu podaci za zadnji element. U tom slučaju neće se više u novi redak 
odštampati zagrade. ' 


Pritiskom na tipku EOF unašaju se upisani podaci u program. Računalo tada analizira is- 
pravnost unesenih podataka i ukoliko ne pronadje nikakvu pogrešku, program se nastavlja slijedećom 
fazom. U slučaju neke pogreške, računalo će ispisati vrstu pogreške i neće nastaviti rad dok se 
pogreška ne ispravi. Kod ovih podataka mogu se pojaviti samo tri vrste pogrešaka i to: 


THE ABOVE SPECIFICATION CONTAINS nae BLOCK NUMBER 
€ 3 ( ) ( 
(Gornja specifikacija sadrži nepropisni broj elemenata) 


: Ova će se poruka odštampati u slučaju kad se broj elementa ne nalazi izmedju 1 i 75. Re- 
dak koji sadrži krivi broj odštampat će se prije ove poruke, a ispod njega odštampat će se četiri 
para zagrada u koje treba unijeti ispravne brojeve pomoću tastature. 


IMPROPER PARAMETER SPECIFICATION FOR ELEMENT 
(Specificirani su parametri koji ne odgovaraju elementu) 

Ova poruka označuje da specificirani parametri nisu u skladu s ranije definiranim elemen- 
tom. Takva pogreška nastaje kad se na primjer definiraju parametri za element koji po svojim oso- 


binama nema pridruženih parametara, ili kad se npr. upiše iznos parametra u kolonu za parametar 3, 
za element koji ima samo jedan ili dva parametra. 


NO CORRESPONDING CONFIGURATION STATEMENT 
(Nema odgovarajućeg elementa u konfiguraciji) 


Ova poruka kazuje da u podacima o konfiguraciji ne postoji element s brojem koji je nave- 
den u podacima o parametrima i početnim uvjetima. 











Niti za vrijeme ove faze nije potrebno ništa s računalom raditi ako su podaci na bušenim 
karticama bez pogrešaka, pa se program nastavlja na slijedeću fazu. 


4.5. FAZA PRIMANJA PODATAKA O GENERATORIMA FUNKCIJA 


Ova faza počinje štampanjem naslova: 


FUNCTION GENERATOR SPECIFICATIONS 
(Podaci za generatore funkcija) 


samo u slučaju ako u podacima o konfiguraciji postoji barem jedan element F (generator funkcije). 
U protivnom slučaju cijela ova faza bit će preskočena i računalo će od faze primanja podataka o 
parametrima automatski prijeći na fazu računanja. 


Nakon štampanja naslova očitavaju se kartice s podacima o generatorima funkcija i ispisu- 
ju se na listu. Podaci se neće ispisati na listu ako je bila uključena preklopka 11. 


Na kraju i ove grupe kartica mora se nalaziti jedna prazna kartica koja će označiti kraj 
ulaska podataka o generatorima funkcija. 


Žele li se unijeti podaci za još koji generator funkcije pomoću tastature, potrebno je 
uključiti preklopku 3 prije nego čitač naidje na praznu karticu. To znači da preklopku 3 treba uk- 
ljučiti najkasnije za vrijeme štampanja posljednjeg retka s podacima za prethodni generator funk- 
cije. Računalo će tada odštampati pet pari zagrada u koje treba upisati podatke kao u prvi red li- 
Ste FUNCTION GENERATOR DATA. Nakon unošenja ovih podataka, računalo odštampa četiri para zagrada 
za drugi red podataka i na kraju tri para zagrada za treći red podataka: 


(ah ( ) ( ) ( 
( ) ( ) ( ) ( 2. 
( ) ( ) ( 


Ukoliko je to posljednji generator funkcija (najviše ih može biti tri), onda je potrebno 
prije unošenja posljednjeg podatka pomoću tastature isključiti preklopku 3, kako računalo ne bi u 
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očekivanju podataka za još jedan generator funkcije počelo ponovno štampati zagrade. 


U slučaju da parametri Pi, i P2 koji označavaju maksimalnu i minimalnu vrijednost ulazne 
varijable nisu bili ranije definirani u listi podataka o parametrima, računalo će zatražiti ove 
podatke ispisivanjem poruke 


SPECIFY UPPER AND LOIHER LIMITS FOR F.G. BLOCK NIH 
(Navedi gornju i donju granicu za generatc» funkeija, element NN) 


Nakon upisivanja numeričkih vrijednosti ovih parametara i uz pretpostavku da je preklop- 
ka 3 bila na vrijeme isključena, program će se pritiskom na tipku EOF nastaviti na slijedeću fazu, 
jasno uz uvjet da su svi podaci bili u redu. 


: Istu će poruku računalo odštampati i u slučaju unošenja podataka o generatorima funkcija 
pomoću bušenih kartica, ako parametri nisu bili u redu uneseni s grupom podataka o početnim uvje- 
tima i parametrima za generatore funkcija. To će biti na primjer slučaj kad parametri zamijene 
svoja mjesta. Naime, parametar 1 (gornja granica) mora biti veći od parametra 2 (donja granica). 


Ukoliko generator funkcija čiji su podaci uneseni u ovoj fazi, nije bio specificiran u 
listi s podacima o konfiguraciji, odštampat će se poruka 


BLOCK NN HAS NOT BEEN SPECIFIED AS A FUNCTIOH GENERATOR 
(Element NN nije bio specificiran kao generator funkcije) 


Nakon toga odštampat će se prvi redak s pet pari zagrada u koje treba upisati novi broj 
elementa koji odgovara broju elementa označenog kao generator funkcije u podacima o konfiguraciji. 
U preostale četiri zagrade treba upisati prva četiri podatka za generator funkcije. Pogreška će se 
moći na ovaj način ispraviti samo u slučaju kad je već generator funkcija bio definiran u podacima 
o konfiguraciji, ali se je zabunon samo u podacima za generatore funkcija upisao krivi broj ele- 
menta. 
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Nakon što su svi podaci u redu primljeni, moguće je u svakom trenutku prekinuti izvršenje 
programa pomoću preklopke 0, i uključenjem preklopke 3 vratiti se u ovu fazu, da bi se prema želji 
promijenio oblik krivulje koja definira pojedini generator funkcije. 


Isto kao i kod prethodne dvije faze, ne će biti potrebno izvršiti nikakve radnje ako su 
podaci o generatorima funkcija ispravno pripremljeni na bušenim karticama. Program će se tada bez 
zaustavljanja nastaviti na slijedeću fazu. 


46. FAZA DAVANJA UPUTA RAČUNALU 


Ovo je prva faza u kojoj nakon pripremnih radova, korisnik mora dati putem tastature upu- 
te računalu u vezi integracije, štampača i crtača krivulje. : 


Kao prvo, računalo pistavlja pitanje u vezi intervala integracije ispisivanjem teksta 


( ) INTEGRATION INTERVAL 
(Interval integracije) 


Bilo gdje unutar ovih zagrada potrebno je pomoću tastature upisati broj koji će predstav- 
ljati interval integracije. Broj mora biti napisan s decimalnom točkom, npr. 5. ili 0.01 itd, Kad 
se kaže unutar zagrada onda se redovno ne piše točno unutar zagrada, nego ustvari redak ispod za- 
grada, ali tako da se broj nalazi na znakovnim mjestima koja su izmedju zagrada. Zagrađe, naime, 
računalo ispisuje na prvom i jedanaestom znakovnom mjestu, a vrijednost za interval integracije 
treba upisati tako da se cijeli broj nalazi na bilo kojim poljima izmedju dva i deset. 


Koji broj odabrati za vrijednost intervala intearacije? Za to nema čvrstog pravila, nego 
treba pravu vrijednost odrediti eksperimentom za svaki problem posebno. Ipak postoje neke smjerni- 
ce. Naime, pristupajući rješavanju nekog zadatka, redovno se vrijednosti rezultata mogu približno 
procijeniti. Istražuje li se na primjer, prijelazni proces nekog regulacijskoa sustava i želi li 
se dobiti oblik i numeričke vrijednosti krivulje vremenskog odziva, onda se za amplitudu i traja- 
nje prijelaznog procesa mogu približno pretpostaviti neke vrijednosti. 
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Za dobivanje vrijednosti intervala integracije, u prvom. redu treba procijeniti trajanje 
prijelaznog procesa. Za prvo računanje neka vrijednost intervala integracije bude 100 do 200 puta 
manja od ukupnog vremena promatranja prijelazne pojave. Na primjer, pretpostavi li se da će pri- 
jelazni proces završiti u roku od 20 jedinica vremena (sekunde, minute, sati, itd.), onda se za 
interval integracije može uzeti npr. vrijednost 0.1. To znači da će računalo provesti 200 prora- 
čuna za traženu krivulju (tj. izračunat će 200 točaka krivulje). 


Uzme li se preveliki interval integracije, integracija ne će biti dovoljno "fina", dok 
će za premali interval integracije cijeli proračun trajati suviše dugo. B 


Drugi kriterij kod odabiranja vrijednosti za interval integracije je očekivana frekvenci- 
ja npr. izlazne veličine. Prema ovoj preporuci, interval integracije treba biti desetak puta manji 
od jedne periode najviše frekvencije. 


Najbolja vrijednost intervala integracije bit će ona za koju će trajanje računanja biti 
najmanje, uz uvjet da točnost ostane u dozvoljenim granicama. Takvu vrijednost moguće je dobiti 
samo eksperimentiranjem. 


Vrijednost intervala integracije mora biti veća od nule. U slučaju da je npr. unutar za- 
grada upisan neki negativni broj ili nula, računalo će odštampati poruku 
INTEGRATION INTERVAL MUST BE GREATER THAN ZERO 


( ) INTEGRATION INTERVAL 
(Interpval integracije mora biti veći od nule) 


Ispod poruke, potrebno je unutar zagrada upisati novu ispravnu vrijednost. Da bi se označio kraj 
unašanja ovog podatka, potrebno je pritisnuti tipku EOF. 


Slijedeće pitanje na koje treba računalu dati odgovor, je u vezi ukupnog vremena. U tu 
svrhu računalo ispisuje 


( ) TOTAL TIME 
(Ukupno vrijeme) 
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Unutar zagrada treba upisati broj koji predstavlja vrijeme, u već ranije definiranim je- 
dinicama (jedinica vremena mora biti ista za sve vremenske veličine u jednom proračunu), za koje 
se želi dobiti proračun, odnosno krivulja neke veličine. Tu je veličinu takodjer potrebno proci- 
jeniti "od oka". Ukoliko je procjena bila loša, računanje se može ponoviti s drugom vrijednošću. 
I ovaj broj mora biti napisan s decimalnom točkom. 


Ukupno vrijeme mora biti veće od intervala integracije. Ukoliko nije, računalo ispisuje 
poruku 


TOTAL TIME MUST BE GREATER THAN INTEGRATION INTERVAL 


( ) TOTAL TIME 
(Ukupno vrijeme mora biti veće od intervala integracije) 


Unutar zagrada potrebno je upisati novu veću vrijednost za ukupno vrijeme. Nakon što je 
podatak upisan, potrebno je pritisnuti tipku EOF, da bi računalo zatražilo podatke u vezi črtača 
krivulje. Kao prvo, računalo ispisuje 


( ) BLOCK FOR Y-AXIS ( ) MINIMUM VALUE ( ) MAXIMUM VALUE 


(Blok za Y-o8, minimalna 1 maksimalna vrijednost) 


Odgovorom na ovo pitanje odredjuje se koja će se veličina grafički iscrtati na crtaču 
krivulja. U tu svrhu potrebno je unutar prve zagrade upisati broj elementa čiji će se izlaz crta- 
ti, a unutar druge i treće zagrade pretpostavljenu minimalnu i maksimalnu vrijednost ove izlazne 
veličine. Ukoliko su minimalna i maksimalna vrijednost krivo pretpostavljene, ili će papir biti 
premali da se na njega nacrta tražena krivulja, ili će krivulja biti suviše malena. 


Prilikom upisivanja treba paziti da je minimalna vrijednost namja od maksimalne. Na prvi 
pogled do zabune ne može doći. Ako se očekuje da će se sve amplitude izlazne veličine nalaziti u 
promatranom intervalu, npr. izmedju 0. i 50., onda je jasno da u drugu zagradu treba upisati 0., 
au treću 50. . Medjutim,očekuju li se amplitude npr. od 0. do -50., onda je minimalna vrijednost 
-50., a maksimalna 0. . 


Ukoliko se ipak upišu krivi podaci unutar zagrada, računalo će odštampati 
ERROR IN Y-AXIS SPECIFICATION 
( ) BLOCK FOR Y-AXIS — ( ) MINIMUM VALUE ( ) MAXIMUM VALUE 


(Pogreška u podacima za Y-0s) 


U zagrade treba tada upisati ispravne vrijednosti. Ista poruka će se odštanpati ako broj 
elementa nije izmedju 1 i 75. 


Nakon unošenja ovih podataka (tj. nakon što se pritisne tipka EOF), računalo traži podat- 
ke za X-os ispisivanjem 


( ) BLOCK FOR X-AXIS ( ) MINIMUM VALUE ( ) MAXIMUM VALUE 
(Blok za X-os, minimalna i maksimalna vrijednost) 

Naime, crtač krivulja (tip 1627) omogučuje dobivanje krivulje u kojoj je jedna varijabla 
funkcija bilo koje druge varijable, tj. dobivanje krivulje tipa y = £(x). U tom slučaju potrebno 
je unutar prve zagrađe upisati broj elementa čija će izlazna veličina predstavljati funkciju X-osi, 
tj. elementa prema kojem će se crtati krivulja y = f(x). U drugu i treću zagradu upisuju se mini- 
malna i maksimalna vrijednost veličine za X-os. I ovdje treba paziti da je maksimalna vrijednost 
veća od minimalne. Ukoliko ipak dodje do zabune prilikom upisivanja, računalo će ispisati 
ERROR IN X-AXIS SPECIFICATION 
( ) BLOCK FOR X-AXIS ( ) MINIMUM VALUE t ) MAXIMUM VALUE 
(Pogreška u podacima za X-os) 


U zagrade treba upisati ispravne vrijednosti. 


Ukoliko se želi dobiti krivulja veličine Y u ovisnosti o vremenu, tj. y = £(t), u tom 
slučaju nije potrebno ništa upisati u zagrade. Računalo će samo automatski postaviti element 76 
(koji predstavlja vrijeme kao nezavisnu varijablu). Minimalna vrijednost bit će nula, a maksimal- 
na će biti jednaka ukupnom vremenu. 


Pritiskom na tipku EOF računalo nastavlja davanje uputa štampanjem teksta 


PREPARE PLOTTER AND PRESS START 


(Pripremi oartač krivulja € pritisni START) 





SET PEN ABOUT ONE INCH FROM RIGHT MARGIN 

(Postavi pisaljku na približno 2,5 em od desnog ruba) 
, Da bi se crtač krivulja tip 
1627 privremio za rad, potrebno je 
postaviti papir preko bubnja koji po- 
miče papir u smijeru X-osi, i postavi- 
ti pisaljku na miesto koje je od de- 
snog ruba papira udaljeno oko 2,5 cm. 
To je potrebno učiniti zato da bi se 
okvir 10x10 inča (25,4x25,4 cm), unu- 
tar kojeg će se crtati krivulja, 
smjestio na sredinu papira koji je 
širok 12 inča (30,5 cm). Promiena va- 
rijable Y pomiče pisaljku horizontal- 
no po nosaču (lijevo - desno). Pi- 
salika, kao i bubanj ne pomiću se 
kontinuirano, nego u koracima po 
1/100 inča (0,254 mm). 


Pritiskom na tipku PROGRAM 
START, crtač krivulja počinje s iscr- 
tavanjem okvira (sl. 13). Crtajući 
Xi Yos, pisaljka na svakih 25,4 mm 
načini križić, tako da na svakoj osi 
načini 10 podjela. Osim toga pisaljka 
još nacrta jednu horizontalnu i jednu 
vrtikalnu crtu po sredini uokvirenog 
dijela papira, tako da se dobiju če- 
tiri jednaka polja 5x5 inča (12,7 x 
12,7 cm). : 


Nakon što je završeno crta- 
nje koordinatnog sustava, koje traje 
cca 40 sekundi, računalo traži sli- 

x jedeći podatak ispisivanjem 





Sl. 13 


( ) PRINT INTERVAL 


(Interval štampanja) 


Unašaniem podatka unutar ovih zagrada odredjuje se koliko će vrijednosti računanih veli- 
čina biti odštampano na štampaču centralne jedinice računala Za prvo računanje bit će dovoljno da 
se odštampa 10 do 20 vrijednosti. Ukoliko se iznosi izlaza iz pojedinih elemenata jako mijenjaju 
i ako treba dobiti njihove iznose u većem broju točaka, potrebno je smanjiti interval štampanja. 


Podatak za interval štampanja dobit će se na temelju ukupnog vremena. Ako je na primjer, 
ukupno vrijeme iznosilo 20, a želi se dobiti 10 odštampanih vrijednosti, interval štampanja mora 
iznositi 2. Za 20 odštampanih vrijednosti, interval štampanja (za ukupno vrijeme 20) mora biti 1, 
za 40 vrijednosti, 0.5 itd. 


Na taj način je interval štampanja u vezi i s intervalom integracije. Prema ranije nave- 
denim preporukama izlazi da je za prvo računanje dovoljno da se odštampa svaka deseta do dvadeseta 
izračunata vrijednost. 


Interval štampanja ne smije biti manji od intervala integracije. U slučaju da je unesena 
manja vrijednost, računalo traži ispravljanje podatka štampanijem 


PRINT INTERVAL CAM NOT BE LESS THAN INTEGRATION INTERVAL 


(Interval štampanja ne može biti manji od intervala integracije) 
( ) PRINT INTERVAL 


Unutar zagrada potrebno ie upisati veću vrijednost intervala štampanja, ili se može sma- 
njiti interval inteqracije..Opet je potrebno pritisnuti tipku EOF da bi računalo prije nego počne 
s računaniem, zatražilo posljednje podatke ispisivanjem 


TIME OUTPUT( ) OUTPUT( ) OUTPUT( ) OUTPUT( ) ONTPUT( ) 
(Vrijeme, izlaz, izlaz, izlaz, izlaz, izlaz) 
Unutar ovih zagrada treba upisati brojeve elemenata za čiju se izlaznu veličinu žele do- 


biti numerički podaci. Istovremeno se kod jednog računanja mogu dobiti podaci za pet različitih 
veličina iz bilo kojeg mjesta u blokovskoj shemi. Ako se u zagradu upiše broj koji se ne nalazi 
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izmedju 1 i 75, ako se upiše broj elementa koji nije specificiran u listi podataka o konfiguraciji, 
ili ako se ne upiše ništa, računalo će u koloni ispod ispisivati nule u obliku 0.0000 . Naime, svi 
izlazni podaci ispisuju se u obliku decimalnog broja koji iza decimalne točke ima četiri decimalha 
mjesta. 


U slučaju kad se u sastavu digitalnoa računala ne nalazi crtač krivulja 1627, program je 
tako sastavljen da se kao rezultat računania, na centralnoj jedinici ispisuju numerički podaci za 
četiri izlaza, a iza četvrtog izlaza štampaju se crtice čiji je broi proporcionalan jednoj od dlz= 
laznih veličina. Na kraju niza crtica odštampa se križić, tako da križići svih redaka prikazuju 
ustvari krivulju, slično kao kad bi crtač krivulja 1627 iscrtavao samo križiće na mjestu koje od- 
govara pojedinom retku numeričkih podataka. 


47. FAZA RAČUNANJA 


U trenutku kad je u posljednju zagradu upisan željeni broj, računalo počinje s računa- 
niem, a crtač krivulja počinje crtati krivulju. Ukoliko se traže numerički podaci za manje od pet 
različitih veličina, onda treba iza unošenja posliednieg broja koji označava broj elementa za či- 
ji se izlaz žele dobiti podaci, pritisnuti tipku EOF da bi računalo započelo proračun. 


Nakon svakoq intervala štampanja ispisuje se jedan redak podataka na štampaču centralne 
jedinice i istovremeno pisaljka crtača krivulja načini na krivulji mali križić, da bi se štampani 
numerički podaci mogli lakše usporediti s nacrtanom krivuliom. . 


Program za crtanje krivulje sastavljen je tako da onemogućava crtanje krivulje izvan 
okvira, iako amplituda varijable na Y-osi prijedje pretpostavljeni iznos. 





4.8. FAZA ZAHTJEVA PO ŽELJI 


Iznad tastature na centralnoj jedinici nalazi se 16 preklopki pomoću kojih se može djelo- 
vati na odvijanje programa. Da se ne bi moralo pamtiti šta se može pomoću pojedine preklopke uči- 
niti, računalo automatski neposredno nakon završetka računanja ispiše za što koja preklopka služi. 


TO SUPPRESS INSTRUCTIONS TURN ON SHITCH 10 
OPTION "SWITCH 

CONF [GURATI ON 

INITIAL CONDITIONS OR PARAMETERS 

FUNCTION GENERATOR INTERCEPTS 

INTEGRATION SPECIFICATIONS 

PRINT INTERVAL 

PRINT VARIABLES 

NEW PLOT FRAME 

PLOTTER X-AXIS * 

PLOTTER Y-AXIS 

SUPPRESS TYPING OF CARD INPUT DATA ll 


DL FUuRNH 


PUNCH UPDATED DATA DECK 12 
INTERROGATE BLOCK OUTPUTS 13 
SAVE STATUS AT INTERRUPT TIME 14 


RESTART AT PREVIOUS INTERRUPT POINT 15 
AFTER SELECTING DESIRED OPTION PRESS START 


Gornje instrukcije računalo će odštampati u svim slučajevima kad se računanje završi ili 
prekine, a to se može dogoditi zbog slijedećih razloga: 


a) ako je dostignuto ukupno vrijeme 
b) ako je uključena preklopka br, 0 iznad tastature na centralnoj jedinici 
c) ako je računanje završeno pomoću elementa Q 


d) uslijed smetnje u procesu računanja. 
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Ako su svi podaci uredno uneseni u računalo, računanje će se obaviti u vremenskom inter- 
valu koji je dan u obliku ukupnog vremena. Nakon tog intervala računanje je gotovo i računalo se 
zaustavlja. 


Ukoliko se računanje prekine prije isteka ukupnog vremena, računalo će ispisati poruku s 
razlogom zašto je došlo do prekida. Ako se želi računanje prekinuti, da bi se npr. izvršile neke 
promjene prema unesenim podacima, potrebno je uključiti preklopku 0. U tom slučaju, kao razlog 
prekida, računalo ispisuje: 


RUN TERMINATED BY SWITCH 0 
(Računanje prekinuto pomoću preklopke 0) 


Nakon što se ispiše ova poruka, potrebno je isključiti preklopku 0. 


Ako je simulaciona shema sastavljena tako da sadrži element Q, računanje će se završiti 
kad prvi ulaz postane veći od drugog. U tom slučaju računalo će odštampati poruku 


RUN TERMINATED BY QUIT ELEMENT 
(Računanje završeno pomoću elementa "gotovo") 


U slučaju kad dodje do poteškoća u procesu računanja, računalo ispisuje 


ERROR IN PROCESSING 


(Smetnja u procesu računanja) 


To se može dogoditi u slučaju kad npr. nisu specificirane gornja i donja granica za gene- 
rator funkcije ili kad npr. drugi ulaz u element za dijeljenje postane jednak nuli. Računanje se 
može nastaviti tek kad se pogreška ispravi. 


Kao što je ranije spomenuto, računalo će nakon jedne od ovih poruka ispisati sve upute u 
vezi preklopki na centralnoj jedinici. Ukoliko se ne želi da računalo ispiše ove upute, potrebno 
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je prije uključivanja preklopke 0, uključiti preklopku 10. Preklopka 10 može se uključiti i kad je 
već počelo ispisivanje uputa. U tom slučaju računale će, nakon što ispiše redak do kraja, prestati 
s daljnjim ispisivanjem. Uključivanjem preklopke 10 ne Će se ispisivati instrukcije u vezi preklo- 
pki, ali će se ispisati samo zadnji redak 


AFTER SELECTING DESIRED OPTION PRESS START 
(Nakon odabiranja željenih preklopki pritisni tipku PROGRAM START) 


Prilikom ispisivanja ovih uputa, računalo će ispisati samo one upute koje su potrebne za 
odvijanje programa koji se upravo računa. Tako će se npr. poruka u vezi preklopke broj 3 odštampa- 
ti samo onda kad simulirani sustav sadrži barem jedan generator funkcije, dok će se poruke u vezi 


preklopki 7, 8 i 9 ispisati samo u slučaju kad je kao izlazna jedinica uključen i crtač krivulja 
1627. - : 


KAKO SE POMOĆU POJEDINE PREKLOPKE MOŽE UTJECATI NA ODVIJANJE RAČUNANJA? 


Preklopka 1 
CONFIGURATION (Konfiguracija) 


Pomoću ove preklopke moguće je dodati, promijeniti ili izbrisati podatke u vezi podataka 
o konfiguraciji. Ukoliko u konfiguraciji postoji i generator funkcije, potrebno je uključiti i 


preklopku 3. Preklopku 1 treba isključiti prije nego se unesu podaci o konfiguraciji za posljednji 
element. B : 


Preklopka_2 


INITIAL CONDITIONS AND PARAMETERS (Početni uvjeti + parametri) 


Uključivanjem ove preklopke moguće je dodati, promijeniti ili izbrisati podatke u vezi po- 
četnih uvjeta i parametara. Preklopku 2 treba isključiti prije nego li se unesu početni uvjeti ili 
parametri za posljednji element. 


Preklopka 3 
FUNCTION GENERATOR INTERCEPTS (Zahvat u vezi generatora funkcija) 
Uključivanjem ove preklopke moguće je dodati, promijeniti ili izbrisati podatke u vezi 


generatora funkcija. Preklopku 3 treba isključiti prije nego što se unesu podaci za posljednji ge- 
nerator funkcije. 


Preklopka_4 
INTEGRATION SPECIFICATIONS (Podaci o integraciji) 


Uključivanjem ove preklopke mogu se promijeniti interval integracije i ukupno vrijeme. 
Postupak je završen kad su podaci uredno uneseni. 


.Preklopka b) 
PRINT INTERVAL  (Interval štampanja) 


Uključivanjem ove preklopke može se promijeniti ranije unesena vrijednost za interval 
štampanja. Postupak je završen kad je novi podatak uredno unesen. 


Preklopka 6 
PRINT VARIABLES (Varijable čiji se iznosi štampaju) 


Uključivanjem ove preklopke mogu se promijeniti brojevi elemenata za čije se izlaze traži 
traži ispisivanje numeričkih vrijednosti. Postupak je završen ka su podaci uredno uneseni. 


Preklopka 7 


NEW PLOT FRAME (Novi okvir za krivulju) 
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Uključivanjem ove preklopke, crtač krivulja 1627 nacrtat će novi okvir za krivulju koju 
treba nacrtati na razmaku od 12,5 cm od prethodnog okvira. Korisnik pri tome ne treba ništa više 
raditi, jer će se okvir automatski smjestiti na sredinu papira, a pisaljka će se sama postaviti u 
početni položaj. š 


Preklopka_8 


PLOTTER X-AXIS (Podaci za *os) 


Uključivanjem ove preklopke može se promijeniti element čiji će izlaz biti predstavljen 
s horizontalnom X-osi, kao i iznos minimalne i maksimalne vrijednosti njegove veličine. Postupak 
je završen kad su svi podaci uredno uneseni. 


Preklopka 9 


PLOTTER Y-AXSIS  (Podaet za Y-os) 


Uključivanjem ove preklopke može se promijeniti element čiji će izlaz biti predstavljen 
s vertikalnom Y-osi, kao i iznos minimalne i maksimalne vrijednosti njegove veličine. Postupak je 
završen kad su svi podaci uredno uneseni. 


Preklopka 10 


TO SUPPRESS INSTRUCTION TURN ON SWITCH 10 (Da se spriječi štampanje uputa uključi pre- 
klopku 10) 


Ukoliko se ne želi da računalo ispisuje upute u vezi preklopki na centralnoj jedinici ra- 
čunala, treba uključiti preklopku broj 10. Ako se predvidja da će tijekom računanja biti više pre- 
kida, npr. zbog promjene ili unošenja novih podataka, onda se preklopka 10 može ostaviti uključena 
za cijelo vrijeme računanja. 





Preklopka 11 
SUPPRESS TYPING OF CARD INPUT DATA (Sprečavanje štampanja ulaznih podataka) 


Ako se uključi ova preklopka, štampač na centralnoj jedinici neće ispisivati ulazne po- 
datke koje čitač kartica očitava s bušenih kartica. U tom slučaju čitanje bušenih kartica ide mno- 
go brže. Ovi se podaci neće štampati čim se uključi preklopka 11, tj. nije potrebno prije nje uk- 
ljučiti preklopku 0. U vezi ove instrukcije korisnik ne mora ništa više raditi. Preklopka 11 može 
se isključiti u bilo kojem trenutku. a 


Preklopka 12 


PUNCH UPDATED DATA DECK (Bušenje kartica prema unesenim podacima) 


Često je puta potrebno tijekom računanja provesti mnoge izmjene kako u konfiguraciji, ta- 
ko i u parametrima i početnim uvjetima. Da bi se dobio paket bušenih kartica prema posljednjim iz- 
mjenama u podacima treba uključiti preklopku 12. Računalo će tada ispisati poruku 


LOAD BLANK CARDS IN CARD READ/PUNCH AMD PUSH START 
(Stavt prazne kartice u čitač/bušač kartica i pritisni tipku PROGRAM START) 

Računalo sada čeka da se stave prazne kartice u čitač/bušač kartica 1442. Nakon toga tre- 
ba pritisnuti tipku PROGRAM START na čitaču kartica 1442, a zatim tipku PROGRAM START na central- 
noj jedinici. Počinje bušenje kartica koje sigurno predstavljaju stanje u trenutku kad su izbuše- 


ne. Izbušit će se kartice s podacima o konfiguraciji, početnim uvjetima, parametrima i generatori- 
ma funkcija. Kad su izbušene sve kartice, računalo se zaustavlja. 


Preklopka 13 


INTERROGATE BLOCK OUTPUTS (Ispitivanje izlaza iz elementa) 


Uključivanjem ove preklopke 13 moguće je dobiti numerički podatak izlaza iz bilo kojeg 
elementa u trenutku prekida. Računalo ispisuje 


OUTPUT INTERROGATION OPTION 
( ) BLOCK 
(Ispitivanje izlaza) 


Unutar zagrada treba upisati broj elementa za čiji se izlaz traži numerička vrijednost. 
Nakon toga računalo ispisuje, npr. : 


OUTPUT OF BLOCK 2 -4.1597 

( ) BLOCK 

(Izlaz iz elementa 2) 

ako je npr. zatražen podatak za element označen brojem 2. Neposredno iza odgovora na postavljeno 
pitanje, računalo ispiše novu zagrađu u koju se može upisati broj slijedećeg elementa za koji se 
želi dobiti vrijednost izlazne veličine. To se ponavlja tako dugo dok je preklopka 13 uključena. 


Postupak je završen kad se isključi preklopka 13 ili kad se pritisne tipka EOF bez da se upiše 
broj unutar zagrada. : 


Ukoliko se unutar zagrada upiše negativan broj ili broj veći od 75, računalo će ispisati 


ENTRY ERROR 
( ) BLOCK 


(Ulazna pogreška) 


Unutar zagrada treba upisati ispravni broj elementa. 


Preklopka 14 


SAVE STATUS AT INTERRUPT TIME (Zadržavanje stanja u trenutku prekida) 


Pomoću ove preklopke moguće je sačuvati stanje u trenutku prekida. Izlazi iz svih eleme- 
nata spreme se na disk, da bi se nakon nekog vremena (za koje se može dalje odvijati računanje) 
moglo ponovno krenuti od mjesta prekida u proračunu. To je posebno interesantno npr. kod ispiti- 
vanja osjetljivosti rješenja na promjene pojedinih parametara, jer u tom slučaju nije potrebno po- 
navljati proračun od početka. i 
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Na disk se može spremiti samo jedno stanje. To znači, ako se preklopka 14 uključi ponovno 
u toku računanja, izbrisat će se prvobitno zapamćeno stanje i upisat će se stanje prilikom drugog 
uključivanja preklopke 14. 


Preklopka 15 
RESTART AT PREVIOUS INTERRUPT POINT (Ponovno kretanje od točke prekida) 


Uključivanjem preklopke 15 program se vraća na situaciju koja je spremljena na disku uk- 
ljučivanjem preklopke 14, (a ne od bilo koje točke prekida kako bi se moglo razabrati iz samog na- 
slova). 


Postupak kod rada s preklopkama 


Ako se želi upotrebljavati preklopke na centralnoj jedinici, potrebno je najprije u ođa- 
branom trenutku zaustaviti odvijanje programa uključivanjem preklopke 0. Računanje se može zausta- 
viti u bilo kojem trenutku. Računalo će ispisati razlog zbog kojeg je došlo do prekida, tj. 

RUN TERMINATED BY SWITCH 0 3 


i odmah nakon toga (ako je bila uključena preklopka 10) 
AFTER SELECTING DESIRED OPTION PRESS START 


Ukoliko preklopka 10 nije bila uključena, odštampat će se izmedju prve i druge poruke cijela lista 
uputa čemu koja preklopka služi, 


Sada se željene preklopke mogu uključiti. Istovremeno se može uključiti više preklopki. 
Bez obzira kojim su redom preklopke uključivane, računalo će izvršavati promjene ovim redoslijedom 


13, 14,1,2,3,12,4,7,9,8,5,6,15. 


Ove se preklopke mogu uključiti i prije nego što je program zaustavljen uključivanjem pre- 
klopke 0. 
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Za vrijeme dok program traži unošenje podataka gori sijalica KB SELECT. 


Nakon što se uključe odabrane preklopke treba pritisnuti tipku PROGRAM START. Računalo 
tada najprije ispiše koje su preklopke bile uključene 


SWITCHES SET ON WERE NN NN 
(Bile su uključene preklopke NN NN... 44) 
gdje umjesto NN stoje brojevi svih uključenih preklopki (osim preklopke 10). 


Tek nakon toga računalo ispisuje tekst koji odgovara radnji pojedinih preklopki. Nakon 
upisivanja potrebnih podataka računalo automatski nastavlja s računanjem. U nekim je slučajevima 
potrebno pritisnuti tipku EOF da bi računalo nastavilo s radom. 


Ako se prilikom upisivanja podataka pomoću tastature načini pogreška (udari krivi broj, 
broj upiše izvan zagrađa ili sl.), potrebno je pritisnuti tipku ERASE FIELD (iabriši polje) i po- 
novno ispisati cijeli redak. Kad se ispiše cijeli redak potrebno je pritisnuti tipku EOF da bi se 
upisani podaci unijeli u program. 


U slučaju kad se nakon zaustavljanja računanja pomoću preklopke 0, ne uključi ni jedna 
preklopka, računalo nakon što se pritisne tipka PROGRAM START, ispisuje 


SMITCHES SET OH WERE * + o* 


U pravilu preklopke treba isključiti odmah nakon što je obavljeno ono što se je pomoću 
pojedine preklopke htjelo učiniti. Ukoliko se do kraja računanja ne namjerava više zaustavljati 
program, onda preklopke mogu ostati uključene; to neće imati nikakvog utjecaja na odvijanje pro- 
grama. Medjutim, ako se program ipak prekine, a npr. nekoliko preklopki je ostalo uključeno, onda 
će računalo tražiti i podatke koje nismo imali namjeru mijenjati. 


Uključivanjem pojedine preklopke program se automatski vraća na odgovarajuće mjesto u 
programu i od njega kreće dalje. 


5. SPECIJALNI ELEMENTI 


, Kad se prilikom sastavljanja simulacione blokovske sheme ukaže potreba za elementom kojeg 
nema medju 25 osnovnih, potrebno je za svaki takav element sastaviti posebni potprogram prema pra- 
vilima. programiranja u FORTRANU. Takav specijalni element (tj. njegov program u FORTRANU) treba 
smjestiti na odgovarajuće mjesto na disku. Ukoliko su već ranije bili definirani neki drugi speci- 
jalni elementi, treba pogledati koji od njih više nije potreban i na njegovo mjesto upisati novi 
specijalni element. Tako spremljeni specijalni element može se kasnije upotrijebiti pozivanjem po- 
moću odgovarajućeg jezičnog simbola kojeg treba upisati u kolonu u koju se upisuje oznaka elementa 
(TYPE) u listi podataka o konfiguraciji. 


Specijalni elementi čine posebnu grupu od pet elemenata čiji su odgovarajući jezični sim- 
boli brojevi 1, 2, 3, 4i5. Za ove elemente ne postoje gotovi upisani potprogrami kao za sve os- 
tale elemente programa 1130 CSMP. Takav specijalni element može biti npr. krivulja histereze, 
funkcija sinusa, kosinusa, eksponencijalna funkcija i dr. Ovih pet specijalnih programa su privid- 
ni programi da bi se zadovoljio LOADER MONITOR SISTEMA: 





Ovim se potprogramima moraju dati slijedeća: imena: 


SPECIJALNI ELEMENT IME 
ze PBIENE....j, | IME. 
#1 SUBI 
#2 SUB2 
#3 SUB3 
#4 SUB4 


#5 . SUB5 
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Sve varijable povezuju se s ovim potprogramima preko zajedničkog (COMMOM) spremnika, dok 
se u pozivnoj (CALL) naredbi ne pojavljuju nikakvi podaci. 


Na kraju ovog poglavlja prikazani su osnovni zahtjevi za definiranje nekih specijalnih 
elemenata. Sve vrijednosti varijabli i upute o dopuštenom postavljanju preklopki mogu se naći u 
poglavlju PROGRAMMING NOTES AND SWITCH LISTINGS priručnika SYSTEM MANUAL. 


Svi podaci spremaju se u dvije grupe u zajednički spremnik. Jedna se zove INTS i upotreb- 
ljavase za spremanje cjelobrojnih varijabli kao što su na primjer brojevi elemenata. Ova grupa mo- 
ra biti dimenzionirana tako da može primiti 587 članova. Svaki član zahtijeva jednu riječ u cen- 
tralnoj memoriji. Druga grupa zove se REALS i sadrži stvarne vrijednosti varijabli, kao što su pa- 
rametri i izlazne veličine iz pojedinih elemenata. Ova grupa mora biti dimenzionirana tako da može 
primiti 395 članova. Svaki član ove grupe zahtijeva dvije riječi u centralnoj memoriji. 


Vrijednosti varijabli koje mogu biti potrebne, nalaze se u tablici 2. 


VARIJABLA UPUTE 
pi a a č — —————==—===S——— 
I Identificira broj elementa 
Izjednači I s INTS(376) 


J Broj elementa na prvom ulazu elementa I 
J se dobiva pomoću naredbe 
Izjednači MTRX2(1) s INTS(76) 
Dimenzioniraj mjesto MTRX2 kao MTRX2(75) 


[2] 
li 


MTRX2 (I) 


K Broj elementa na drugom ulazu elementa I 
K se dobiva pomoću naredbe . K 
Izjednači MTRX3(1) s INTS(151) 
Dimenzioniraj mjesto MTRX3 kao MTRX3(75) 


MTRX3 (I) 


LI 


L Broj elementa na trećem ulazu elementa I 
L se dobiva pomoću naredbe L 
Izjednači MTRX4(1) s INTS(226) 
Dimenzioniraj mjesto MTRX4 kao MTRX4(75) 


MTRX4 (I) 


_ KA _ 


C(J) Sadrži tekuću vrijednost prvog ulaza elementa I 
C(K) Sadrži tekuću vrijednost drugog ulaza elementa I 
C(L) Sadrži tekuću vrijednost trećeg ulaza elementa I 


Izjednači C(1) s REALS(2) 
Dimenzioniraj mjesto C kao C(76) 


PAR1 (I) Sadrži parametar Pl elementa I 
Izjednači PARI(1) s REALS(81) 
Dimenzioniraj mjesto PARI kao PAR1(75) 


PAR2 (I) Sadrži parametar P2 elementa I 
Izjednači PAR2(1) s REALS(156) 
Dimenzioniraj mjesto PAR2 kao PAR2(75) 
PAR3 (I) Sadrži parametar P3 elementa I 
Izjednači PAR3(1) s REALS(231) 
Dimenzioniraj mjesto PAR3 kao PAR3(75) 
C(I) Spremi izlaz iz elementa I u ovu lokaciju 


€ Nezavisna varijabla vrijeme 
Izjednači T s REALS(77) 


DT B Interval integracije 
Izjednači DT s REALS(78) 


TTOT Ukupno vrijeme proračuna 
Izjednači TTOT s REALS (79) 


TABLICA 2 
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Komplet bušenih kartica naredbi DIMENSION (dimenzioniraj), COMMON (zajednički) i EQUIVA- 
LENCE (zzjednači) sadrži kartice: 


REAL REALS (395) 

INTEGER INTS (587) 

DIMENSION C(76) ,MTRX2(75) ,MTRX3 (75) ,MTRX4 (75) 
DIMENSION PAR1(75),PAR2(75),PAR3(75) 
COMMON  REALS, INTS 

EQUIVALENCE  (REALS (2) , C(1)) 
EQUIVALENCE  (REALS (81), PARI(1)) 
EQUIVALENCE  (REALS(156), PAR2(1)) 
EQUIVALENCE  (REALS(231), PAR3(1)) 
EQUIVALENCE  (INTS(76) , MTRX2(1)) 
EQUIVALENCE  (INTS(151) , MTRX3(1)) 
EQUIVALENCE  (INTS(376) , I ) 


Ispred FORTRAN programa kojim se definira novi specijalni element mora se nalaziti skup 
bušenih kartica koje će poništiti prethodni potprogram na disku. Nakon njih dolazi grupa kartica 
koje omogućavaju spremanje novog potprograma na disk. Teeba napomenuti da kartice DELETE (iabri- 
ši), SUBROUTINE SUBn (potprogram SUBn) i STORE (spremi) moraju imati isto ime kojim je označen do- 
tični potprogram (SUBI, SUB2, SUB3, SUB4 ili SUBS5). 


Dok se ne stekne dovoljno iskustva, treba biti pažljiv kod sastavljanja ovih potprograma. 
Isto tako treba dobro paziti kad se sastavlja program za element koji sadrži memoriju ili spremnik 
dobivenih podataka. U takvom slučaju potrebno je definirati posebna područja za spremanje ovih po- 
dataka unutar potprograma za specijalni element, jer za njih nije predvidjeno uskladištenje u za- 
jednički (COMMON) spremnik. Takav element može biti upotrebljen samo jednom u bilo kojoj konfigu- 
raciji. Nakon što je specijalni elemeut jednom bio upotrebljen za vrijeme računanja, njegovi poda- 
ci će se poništiti. 
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Ako je u proračunu potrebno nekoliko jednakih specijalnih elemenata, tađa treba definira- 
ti toliko specijalnih elemenata koliko je potrebno (ali najviše njih pet), primijenjujući jednake 
potprograme koji se razlikuju jedino u imenu. 


Izuzetak od ovog pravila čini element koji zahtijeva tri ili manje lokacija u spremniku. 
U tom slučaju mnogo je prikladnije upotrijebiti parametre Pl, P2 i P3 za spremanje podataka koji 
pripadaju tom elementu. Tada se pojedini element može upotrijebiti više puta 'u jednoj konfiguraci- 
ji. Na primjer, u potprogramu za generator impulsa s vremenski promjenljivom frekvencijom, podaci 


za ovaj potprogram spremljeni su u parametru 1 bilo kojeg bloka koji upotrebljava ovaj specijalni 
element. : 


Potprogrami nekih specijalnih elemenata 


Potprogram za specijalni element koji sadrži funkciju sinusa 


Er FORTRAN 


Proqramer 
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Potprogram za specijalni element koji sadrži funkciju kosinusa 


[5% FORTRAN 
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Potprogram za specijalni element koji sadrži funkciju arkus tangens 


[57 FORTRAN 
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Potprogram za specijalni element koji sadrži eksponencijalnu funkciju 


L&7 FORTRAN j 


Proqramer Problem Str. 3 





Oznaka Datum Od 
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Potprogram za specijalni element koji sadrži generator impulsa promjenljive frekvencije 
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6. PRIMJERI 


PRIMJER 1 


Mehanički sustav koji se sastoji od mase, opruge i viskoznog odbojnika i na koji djeluje 
sila F, može se opisati diferencijalnom jednadžbom ' 


MX + DX + Kx = K F(t) (1) 


Ako se pretpostavi da se pod mehaničkim sustavom smatra kotač automobila s pripadnim 
amortizerom i oprugom prema slici 14, onda je 


M masa kotača, osovine i drugih pokretnih dijelova 
o konstanta prigušenja amortizera 

K konstanta krutosti opruge 

x pomak iz ravnotežnog položaja 

X brzina 

g ubrzanje 

F(t) sila 





Jednadžba (1) vrijedi uz pretpostavku da je masa 
automobila na koju je učvršćen ovaj sustav, nepomična u ver- 
[rta tikalnom sjeru, ida se opisani sustav može gibati samo u 
SI. 14 jednom smjeru. 


= 82 _ 


Jednadžba (1) može se napisati i u okliku 


žo+ 2Euk + wx = 02 F(t) (2) 

gdje je -D 
; 2čuw M (3) 

2 K 

W= 
M (4) 


i gdje £ predstavlja faktor prigušenja, a w vlastitu frekvenciju sustava. 


Ako se žele istražiti vlastita svojstva sustava masa-opruga-viskozni odbojnik, onda treba 
postaviti da je desna strana jednadžbe jednaka nuli, pa ostaje ' 


MX + DX + Kx = 0 7 - ' (5) 


. To je ujedno poznati oblik diferencijalne jednadžbe tzv. autonomnog sustava drugog reda, 
tj. dinamičkog sustava drugog reda na kojeg ne djeluju vanjske pobude. 


Najjednostavniji način na koji se može promatrati vlastito ponašanje sustava je taj, da 
se masa izvede iz ravnotežnog položaja za npr. pomak x, i naglo pusti. Masa će se vratiti u ravno- 
težno stanje na jedan od tri moguća načina, što ovisi o iznosu i međjusobnom odnosu koeficijenata 
M, DiK. Kik. 


Za vrijednost faktora prigušenja manju od jedan, krivulja odziva imat će oscilatorno pri- 
gušeni oblik; za vrijednosti £ izmedju 0,707 i 1 imat će aperiodsko prigušeni oblik s približava- 
njem s gornje strane, tj. masa će pri vraćanju u početni položaj prijeći ravnotežno stanje, a za- 
tim će se aperiodski približiti konačnom položaju; za £ jednako jedan, krivulja će se aperiodski 
približavati ravnotežnom stanju; za £ veće od jedan, krivulja će se takodjer aperiodski približa- 
vati ravnotežnom stanju ali mnogo sporije. 


Da bi se jednadžba (5) nacrtala u obliku simulacione blokovske sheme sa standardnim ele- 
mentima programa 1130 CSMP, prikladno ju je napisati u takvom obliku da na lijevoj strani ostane 
samo varijabla s najvišom derivacijom. 


ža (DŽ + Xx) (6) 
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Kako sastaviti blokovsku simulacionu shemu? Kao prvo može se uočiti da pored varijable x 
postoji i njezina prva i druga derivacija - dakle moraju se upotrijebiti dva integratora. Ulaz u 
prvi integrator treba biti *. Veličinu X dobit će se ako se izlaz iz elementa pomoću kojeg će se 
zbrojiti Dx i Kx pomnoži s veličinom -1/M. Budući da kod sastavljanja blokovske sheme treba nasto- 
jati da se upotrijebi što manji broj elemenata, treba po mogućnosti sve parametre (ili barem što 
više njih) pridružiti upotrebljenim elementima. Tako je bolje u ovom slučaju umjesto elementa za 
zbrajanje (oznaka +) upotrijebiti ele- 
ment W kod kojeg postoji mogućnost da 
se prije zbrajanja svaki ulaz još i po- 
množi s proizvoljnom konstantom. Umje- 
sti elementa za dijeljenje s parametrom 
M, smještenog iza elementa W, može se 
ući na drugi (ili treći) ulaz integra- 
tora i pomoću parametra P2 pomnožiti 
izlaz iz elementa W s veličinom -I/M. 





Na slici 15 nacrtana je 
shema sa svega tri elementa, dok je 
npr. za shemu na slici 16 (koja je ta- 
kodjer dobra što se tiće rezultata dobivenih na kraju računanja) potrebno čak tri puta više, tj. 
devet elemenata. 


Sl. 15 


Iz ovog primjera 
crtanja blokovske sheme 
vidi se da se pažljivim 
crtanjem može uštedjeti 
veliki broj elemenata, što 
je posebno interesantno 
kod složenih sustava kod 
kojih se može dogoditi da 
se potroši svih 75 raspo- 
loživih elemenata prije 
nego se nacta kompletna 
simulaciona shema. Tada 
treba pokušati smanjiti 
broj elemenata na sličan 
način kao i u ovom primje- 
ru, 





Sl. 16 
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Na osnovi simulacione blokovske sheme (sl. 15) sastavljaju se lista s podacima o konfigu- 
raciji i lista podataka o parametrima i početnim uvjetima. Budući da računalo ne može raditi s op- 
ćim brojevima, potrebno je svaki zadatak zadati s numeričkim vrijednostima. U ovom slučaju neka je 
na primjer zadano: M =2. K = 400. D=14.14 


1130 CONTINUOUS SYSTEM MODELING PROGRAM 


PRIMJER 4 mre 31.72. 


PAGE 


PROGRAM 


ci 


oF 





CODED BY 
CONFIGURATION DATA 


OUTPUT NAME INPUT 4 INPUT 2 INPUT 3 
48 50 





INITIAL CONDITIONS AND PARAMETER DATA 


IC/PAR NAME BLOCK t/c PAR 1 -PAR 2 PAR 3 
ja 6 
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Prema ovim listama treba izbušiti kartice na mjestima na kojima šu upisani podaci. Pored 
toga potrebno je prije samog rada na digitalnom računalu pripremiti još niz podataka koje će tre- 
bati dati računalu za vrijeme samog računanja. 


Kao prvo treba procijeni o trajati prijelazni proces. Pretpostavimo da će 


trajati 2 sekunde. Premda preporuci može se dobiti interval integracije, ako se ukupno vrijeme po- 
dijeli sa sto. Dakle interval integracije neka Bude U. s. z 
—— 


Potrebno je odrediti još jednu vremensku veličinu. To je interval štampanja izlaznih po- 
dataka. AKG smatramo da će >iti dovoljno da se odštampa dvadeset podataka onda iaterval štampanja 
mora iznositi 0.1 s. Budući da se istovremeno mogu dobiti numerički podaci iz pet različitih ele- 
menata, može se kod ovog zadatka zatražiti ispisivanje vrijednosti izlaza na kraju svakog interva- 


la štampanja, za sva tri elementa. 





Ako se želi nacrtati npr. vremenski dijagram puta (tj. krivulja u kojoj je put funkcija 
vremena) kojim se masa vraća u ravnotežni položaj, potrebno je upisati u zagradu bloka za Y-os 
brojku 3, kao i pretpostavljenu minimalnu i maksimalnu vrijednost izlaza iz elementa broj tri. U 
ovom primjeru nije teško zaključiti koje vrijednosti treba upisati. Ako se želi krivulja smjestiti 
na papir tako da ravnotežni položaj bude na sredini nacrtanog okvira, onda treba za minimalnu vri- 
jednost upisati -I., jer je za -1. masa izvedena iz ravnotežnog položaja (što odgovara jediničnoj 
ulaznoj funkciji), dok se za maksimalnu vrijednost treba upisati 1., jer je to maksimalno nadviše- 
nje koje može nastati u slučaju čisto oscilatornog odziva (tj. kač je prigušenje jednako nuli + 
rub stabilnosti). 


Budući da se crta vremenski dijagram, X-os je vrijeme, pa kod traženja podataka za X-os 
ne treba ništa upisati. 


Sada se može pristupiti radu s đigitainim računalom, 


Kao prvo treba provjeriti đa li je stavljen disk broj 5002. Ako treba promijeniti disk, 
potrebno je nakon stavljanja diska pročekati 90 sekundi da disk dobije potrebnu brzinu. Na cen-' 
tralnoj jedinici potrekno je pritisnuti tipke FILE SWITCH, RUN, IMM STOP i RESET, Ove se tipke na- 
laze pored tastature. 


U čitač kartica treba staviti komplet bušenih kartica licem prema dole i to tako da odre- 
zani dio buće s lijeve strane, i pritisnuti tipku START (na čitaču kartica). 
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Komplet bušenih kartica sastoji se u ovom primjeru od slijedećih kartica: 


- kartica COLD START 

- kartica // XEQ CSMP 

- tri kartice s podacima o konfiguraciji 

- jedna prazna kartica 

_ tri kartice s podacima o parametrima i početnim uvjetima 
- barem dvije prazne kartice 


Da u zadatku postoji i generator funkcija, bila bi na kraju dovoljna samo jedna prazna 
kartica. To, dakako, vrijedi samo u slučaju kađ se unaprijed znade da neće biti potrebne nikakve 
izmjene i promjene podataka o konfiguraciji, početnim uvjetima i parametrima i podataka o genera- 
torima funkcija. U protivnom za svaku takvu promjenu potrebna je jedna prazna kartica. 


Ukoliko je stavljeno premalo praznih kartica ua kraju kompleta bušenih kartica, računalo neće do-' 
pustiti ulazak podataka pomoću tastature tako dugo. dok se u čitač kartica ne stavi još praznih 
kartica. 


Na centralnoj jedinici treba pritisnuti tipku PROGRAM LOAD. Svijetle sijalice FILE READY 
i RUN. Računalo očita prve dvije kartice i nakon po prilici 3,5 minute počne odvijanje programa 
ispisivanjem naslova, uputa i ostalog. : 


CONTINUOUS SYSTEM sMODELING PROGRAM 
A DIGITAL ANALOG SIMULATOR "ROGRAM FOR THE IBM 1130 


INSTRUCTIONAL COMMENTS MAY BE SUPPRESSED. AT ANY TIME BY TURNING ON SWITCH 10 

TURN ON SWITCH 1 TO ENTER OR MODIFY CONFIGURATION STATEMENTS VIA THE KEYBOARD 

TURN ON SWITCH 2 TO ENTER OR MODIFY INITIAL CONDITIONS OR ELEMENT PARAMETERS VIA THE KEYBOARD 
TURN ON SWITCH 3 TO ENTER OR MODIFY FUNCTION GENERATOR INTERCEPTS VIA THE KEYBOARD 


CONFIGURATION SPECIFICATION 


OUTPUT NAME BLOCK TYPE INPUT 1 INPUT 2 INPUT 3 
SILA 1 Va 2 3 0 
BRZINA 2 I 9 1 0 


PUT 3 i 2 0 : 0 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PARI PAR2 PAR3 

KOEFICIJENTI 1 14.0400 400.0000 0.0000 

MASA 2 0.0000 -0.5000 0.0000 

POCETNI POMAK 3 -1.0000 0.0000 0.0000 

(0.02 ) INTEGRATION INTERVAL 

( 2. ) TOTAL TIME 

( 3) BLOCK FOR Y-AXIS ( =; ) MINIMUM VALUE ( lo ) MAXIMUM VALUE 
( ) BLOCK FOR X-AXIS ( ) MINIMUM VALUE ( ) MAXIMUM VALUE 


PREPARE PLOTTER AND PRESS START 
SET PEN ABOUT ONE INCH FROM RIGHT MARGIN 


Nakon upisivanja podataka (pomoću tastature) za interval integracije, ukupno vrijeme, za 
Y-os i za X-os, potrebno je bilo svaki put pritisnuti tipku EOF. Nakon upozorenja da se pripremi 
crtač krivulja za rad, treba uključiti crtač krivulja, postaviti papir na valjak i dovesti pisalj- 
ku u položaj cca 2,5 cm od desnog ruba papira, kad se gleda prema papiru. Ovu pripremu crtača kri- 
vulja moglo se je obaviti i ranije, npr. za vrijeme prepisivanja programa iz diska u centralnu me- 
moriju, da se uštedi na vremenu. 


Kad je crtač krivulja pripremljen za rad, treba pritisnuti tipku PROGRAM START, a nakon 
toga unijeti u odštampane zagrade podatak za interval štampanja i brojeve elemenata čije se nume- 
ričke vrijednosti izlaza žele dobiti. Iza unašanja intervala štampanja i brojeva izlaza za Štam- 
panje, bilo je potrebno pritisnuti tipku EOF. 
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Neposredno nakon toga počinje računanje, crtanje krivulje (sl. 17) i ispisivanje podataka. 


« 0.1) PRINT INTERVAL 


TIME OUTPUT( 3) OUTPUT( 2)  QUTPUT(C 1)  OUTPUT( )  QUTPUT( ) 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.9000 
0.100 -0.3039 9.9782 19.5394 0.0000 0.0000 
0.200 0.4055 2.5448 198.1745 0.0000 . 0.0000 
0.300 0.2502 -4.1703 41.1038 0.0000 0.0000 
0:400 -0.1320 -2.2926 -85.2312 0:0000 0.0000 
0.500 -0.1545 1.4295 -41.5890 0.0000 0.0000 
0.600 0.0244 1.4718 30.5583 0.0000 0.0000 
0.700 0:0808 -0.3150 27.8799 0.0000 0.0000 
0.800 0.0088 -0.7912 -7.5487 0.0000 0.0000 
0.900 -0.0368 -0. 0496 -15.4161 0.0000 0.0000 
1.000 -0.0137 0.3695 -0.2374 0.0000 0.0000 
1.100 0.0143 0.1182 7.3851 0.0000 0.0000 
1.200 0.0102 -0.1483 1.9981 0.0000 0.0000 
1.300 -0.0043 -0.0949 -3.0572 0.0000 0.0000 
1.400 -0.0060 0.0474 -1. 7447 0.0000 0.0000 
1.500 0.0005 0.0579 1.0316 0.0000 0.0000 
1.600 0.0031 -0.0081 1.1058 0.0000 0.0000 
1.700 0.0005 -0.0301 -0.2162 0.0000 0.0000 
1.800 -0.0013 -0.0037 -0.5892 0.0000 0.0000 
1.900 -0.0006 0.0136 -0. 0458 0.0000 0.0000 
2.000 0.0005 0.0053 0.2729 0.0000 _ 0.0000 


AFTER SELECTING DESIRED OPTION PRESS START 


Iako je cijeli program sastavljen potpuno korektno i prema svim uputama firme IBM, ipak 
se dobiva utisak, promatrajući nacrtanu krivulju, da nešto nije bilo u redu. Naime, početak krivu- 
lje, kao i vrh prvog nadvišenj: su posve okrugli. Na temelju toga može se izvesti zaključak: ako 
krivulja nije niti po izaledu onatva kakva se očekuje, onda vrlo vjerojatno niti numerički podaci 
nisu ispravni. Mečjutim, ako se zu nešto kaže da nije točno, onda je prije toga bilo potrebno de- 
finirati šta se smatra pod "dovoljno točno", tj. treba postaviti granice unutar kojih se smatra da 
rezultati zadovoljavaju. 









: Jedino "loše" mjesto u pro- 
vedenom proračunu je integracija. 
Buđući da postupak integracije nije 
kontinuirani nego se provodi u kora- 
cima, može se točnost računanja po- 
visiti smanjenjem koraka integraci- 
je, tj. smanjenjem intervala integra- 
cije. U tu svrhu provedeno je još 
nekoliko proračuna kod kojih su svi 
podaci osim intervala integracije 
ostali isti. 


Zbog potpunosti proučavanja 
kako utječe interval integracije na 
točnost dobivenih rezultata, prove- 
dena su najprije tri proračuna s još 
većim intervalima integracije. Narav- 
no da su i dobiveni rezultati lošiji. 


Kod ovakvog eksperimentira- 
nja moguće je uključivanjem odredje- 
nih preklopki iznad tastature na cen- 
tralnoj jedinici, vratiti se na odre- 
djeni dio programa, promijeniti neke 
podatke i provesti ponovni proračun 
do kraja zajedno s crtanjem nove kri- 
vulje. 


Prema gore postavljenom zah- 
tjevu treba promijeniti interval in- 
tegracije, nacrtati novu krivulju 
(na novi dio papira) i dobiti nove 
numeričke podatke u vezi izlaza iz 
pojedinih elemenata blokovske sheme. 


Računalo je s time u vezi 
već ispisalo poruku nakon završenog 
prvog proračuna: AFTER SELECTING 


1 


Si. 18 





DESIRED OPTION PRESS START). U tu 
svrhu porebno je uključiti: 


-  preklopku 4 - da se omogući una- 
šanje novog intervala integracije 


- preklopku 7.- da se nacrta novi 
okvir za sliku. (Ukoliko se želi 
nacrtati novu krivulju na pret- 
hodnu sliku, ovu preklopku nije 
potrebno uključivati. Medjutim 
zbog lakše analize dobivenih kri- 
vulja zatraženo je kod svih ovih 
primjera da se svaki put nacrta 
nova slika.) 


-  preklopku 10 - da se spriječi ne- 
potrebno ispisivanje raznih uputa 
više puta. (Ovu se preklopku ta- 
kodjer nije moralo uključiti, ali 
se na taj način štedi vrijeme.) 


Nakon što su uključene pre- 
klopke 4, 7 i 10, treba pritisnuti 
tipku PROGRAM ŠTART. Kao prvo, raču- 
nalo najprije ispiše koje su preklop- 
ke bile uključene, a zatim zatraži 
podatak za interval integracije (ko- 
ji će u ovom slučaju biti samo 50 
puta manji od ukupnog vremena, tj. 
0.04), zatim podatak za ukupno vri- 
jeme. Ukupno vrijeme će ostati za. 
sve primjere, u kojima se proučava 
utjecaj različitih intervala inte- 
gracije, isto, tj. jednako 2. . 


Iza toga crta se novi okvir 
sa sliku, pisaljka se sama postavlja 
1 ishodište koordinatnog sustava i 
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program se dalje odvija automatski do kraja proračuna bez ikakve intervencije operatera. 


SWITCHES SET ON WERE ko7 
( 0.04 ) INTEGRATION INTERVAL 
( 2. ) TOTAL TIME 


TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 QUTPUT_ 0 


0.000 “-1.0000 0.0000 -400.0000 0.0000 0.0000 
0.120 -0.0555 9.1216 106.774 0.0000 0.0000 
0.200 0.4242 1.6492 192.9825 0.0000 0.0000 
0.320 0.0988 -4.3092 -21.4266 0.0000 0.0000 
0.400 -0.1663 -1.4952 -87.6656 0.0000 0.0000 
0.520 -0.0774 1.9161 -3.8765 0.0000 0.0000 
0.600 0.0582 0,9957 37.3628 0.0000 0.0000 
0.720 0.0486 -0. 7967 8.1905 0.0000 0.0000 
0.800 -0.0165 -0,5764 -14.7417 0.0000 0.0000 
0.920 -0.0272 0.3042 -5.4796 0.0000 0.0000 
1.000 0.0022 0.3053 5.2111 0.0000 0.0000 
1.120 0.0140 -0.1019 4.1777 0.9000 0.0000 
1.200 0.0016 -0.1510 -1.5073 0.0000 0.0000 
1,320 -0.0068 0.0260 -2.3516 0.0000 0.0000 
1.400 -0.0019 0.0702 0.2292 0.0000 0.0000 
1.520 0.0031 * -0.0013 1.2203 0.0000 0.0000 
1.600 0.0014 -0.0307 0.1211 0.0000 0.0000 
1.720 -0.0013 -0.0045 =0.5931 0.0000 0.0000 
1.800 -0.0008 . 0.0125 -0.1598 0.0000 0.0000 
1.920 0.0005 0.0043 0.2714 0.0000 0.0000 


2.000 0.0005 -0.0047 0.1185 0.0000 0.0000 





Iz slike 18 vidi se da je 
situacija lošija od prethodne, što 
je bilo i za očekivati. Koliko je 
lošija može se zaključiti mjerenjem 
duljine odsječaka na vremenskoj osi 
gdje krivulja siječe vremensku os, 
kao i usporedjivanjem iznosa poje- 
dinih maksimuma krivulje s odgova- 
rajućim iznosima krivulje na slici 
17. 


Još bolja usporedba može se 
provesti ako se pogledaju iznosi 
krivulja usakon svakih 0.1 s, što je 
označeno križićima na samoj krivulji, 
a čemu odgovaraju i numeričke vrijed- 
nosti odštampane u drugoj koloni (= 
izlaz iz elementa broj 3). 


Ako se (kao u ovom primjeru) 
želi nekoliko puta promijeniti istu 
veličinu potrebnu za proračun, onda 
su preklopke 4, 7 i 10 mogle ostati 
cijelo vrijeme uključene. Da bi se 
dobila slijedeća krivulja i slijede-. 
ći proračun, potrebno je nakon zavr- 
šetka prethodnog proračuna samo pri- 
tisnuti tipku PROGRAM START i unije- 
ti novu vrijednost intervala inte- 
gracije. 


Budući da se promjena in- 
tervala integracije traži s istom 
preklopkom kao i promjena ukupnog 
vremena, potrebno je svaki put po- 
novno upisati i iznos ukupnog vre- 
mena. 
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Za slijedeći primjer odabran je još veći interval integracije (koji je samo 20 puta manji 
od ukupnog vremena) i koji je u ovom slučaju upravo jednak intervalu štampanja. To znači, da će 
svaka vrijednost koju računalo izračuna biti i odštampana (sl. 19). 


) INTEGRASION INTERVAL 


m 


( 2. ) TOTAL TIME 


TIME OUTPUT _ 3 OUTPUT _2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.0000 12.9300 182.8301 0.0000 0.0000 
0.200 0.8359 -5.9100 250.8027 0.0000 0.0000 
0.300 -0.3821 -8.1072 -267.4677 0.0000 0:0000 
0.400 -0.5241._ < 8.6459 -87.3992 0.0000 0.0000 
0.500 0.5590 2.8252 263.5308 0.0000 0.0000 
0.600 0.1826 -8.5186 -47.3942 0.0000 0.0000 
0.700 -0.5507 1.4320 -198.6291 0.0000 0.0000 
0.800 0.0990 6.4207 130.4063 0.0000 0.0000 
0.900 0.4151 -4.2154 106.4334 0.0000 0.0000 
1.000 -0.2725 -3.4405 -157.6578 0.0000 0.0000 
1.100 -0.2224 5.0963 -16.9089 0.0000 0.0000 
1.200 0.3295 0.5466 139,5186 0.0000 0.0000 
1.300 0.0353 -4.5099 -49.6359 0.0000 0.0000 
1.400 -0.2916 1.6045 -93,9396 0.0000 0.0000 
1.500 0.1037 3.0366 84.4293 0.0000 0.0000 
1.600 0.1963 -2. 7292 39.9359 0.0000 0.0000 
1.700 -0.176L -1.2909 -88.8302 0.0000 0.0000 
1.800 -0.0835 2.8714 - 7.2186 0.0000 0.0000 
1.900 0.1856 -0.2333 70.9558 0.0000 0.0000 


2.000 —040151 -2.2936 -38.4663 0.0000 0.0000 
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Slijedeći primjer (sl. 20) pokazuje do kakvih besmislica može dovesti daljnje povećavanje 
intervala integracije, koji u ovom slučaju iznosi samo deseti dio ukupnog vremena i koji je veći 
čak i od intervala štampanja numeričkih rezultata. 


( 0.2) INTEGRATION INTERVAL 


( 2. ) TOTAL TIME 


TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.09000 0.0000 -400.0000 0.0000 0.0000 
0.200 3.C090 11.7200 1365.7204 0.0000 0.0000 
0.400 -8.5152 “15. 1757 -4388.2627 0.0000 0.0000 
0.600 20.5543 354. 1030 13220.7304 0.0000 0.0000 
0.800 -40,6524 -1449,67G62  -36839.535750 0.0000 0.0000 
1.000 37.C061 5428.4189 91800.2345 0.0000 0.0000 
1.200 205.2875  -18975.0000 -186191.4066 . 0.0000 0.0000 
1.400 -1727.7980 62380.3907  190939.3754 0.0000 0.0000 
1.600 8838,6847 -192735.6254_ 810272.2514 0.0000 0.0000 
1.800  -37810.9532  554455.5264-7284238.0117 0.0000 0.0000 
2,000  145924,5003«“1449967,752937867256.1093 0.0000 0.0000 


Ako se pogledaju numerički podaci, onda se vidi da se je računalo "snašlo" i ispisalo sa- 
mo one podatke koje je izračunalo. Naime, iako je zatraženo da se ispisuju rezultati nakon svake 
0.1 s, računalo je ispisalo rezultat svakih 0.2 s, jer na polovici intervala integracije, numerič- 
ke vrijednosti izlaza iz pojedinih elemenata nisu definirane. 


Sada se može provesti nekoliko proračuna s manjim intervalom integracije nego što je bio 
u prvom slučaju, pa da se pokuša zaključiti do kuda ima smisla ići sa smanjivanjem tog intervala. 





Kumi 
Lii 


Primjer izračunan s intervalom integracije 0.01, #3. intervalom koji je 200 puta manji od 
ukupnog vremena prikazuju slika 21. i siijedeći rumerički rezvitati: 


( 0.01 ) INTEGRATION INTERVAL 


( 2. ) TOTAL TIME 


TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0600 -400.0000 0.0000 0.0000 
0.100 -0.3147 10.02753 15.9015 0.9000 0.0000 
0.200 0.4037 2.753 200.+609 0.0000 0.0000 . 
0.300 0.2653 -4.1575 47.3229 0.0000 0.0000 
0.400 -0.1250 -2.4947 -55.2586 0.0000 0.0000 
0.500 -0. 1644 1.3522 -46.6412 0.0000 0.0000 
0.600 0.0161 1.4948 28.9713 0.0000 0.0000 
0.700 0.0850 -0.224U4 30.8304 0.9000 0.0000 
0.800 0.0155 -0.8435 =5,1253 0.9000 0.0000 
0.900 -0.0374 -0.1219 -16.6880 0.0000 0.0000 
1.000 -0.0179 0.3800 -1.7818 0.0000 0.0000 
1.100 0.0134 0.1642 7.6909 0.0000 0.0000 
1.200 0.0125 -0.0411 2.9881 0.0000 0.0000 
1.300 -0.0032 -0.1193 -2.9486 0.0000 0.0000 
1.400 -0.0070 0.0364 -2.2756 0.9000 0.0000 
1.500 -0.0004 0.0685 0.8198 0.0000 0.0000 
1.600 0.0033 0.0010 1.3410 0.0000 0.0000 
1.700 0.0011 0.0333 -0.0332 0.0000 0.0000 
1.800 -0.00153 -0.0096 -0.6665 0.9000 0.0000 
1.900 -0.0009 0.0137 -0.1670 0.0000 0.0000 
2.000 0.0004 0.0085 0.2808 0.0000 0.0000 


Krivulja dobivena na slici 21 po obliku izgleda dovoljno dobra. Usporedjujući je s krivu- 
ljom na slici 17 (interval integracije 0.02), može se vidjeti da se vrijednosti razlikuju za cca 
2%. 





= 0 = 


Smanji li se interval integracije na petstoti dio ukupnog vremena (sl. 22), poboljšanje 
prema prethodnom slučaju iznosi samo cca 1%. 


( 0.004 ) INTEGRATION INTERVAL 


( 2, ) TOTAL TIME 
TIME OUTPUT 5 OUTPUT 2 OUTPUT 1 OUTPUT _ 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.3176 10.0452 14.9954 0.0000 0.0000 
0.200 0.4037 2.8139 201.2509 0.0000 0.0000 
0.300 0.2695 4.1603 48.9876 0.0000 0.0000 
0.400 -0.1233 “2.5516 -85.4189 0.0000 0.0000 
0.500 -0.1673 1.3347 -48.0603 0.0000 0.0000 
0.600 0.0139 1.6309 28.6168 0.0000 0.0000 
0.700 0.0863 -0.2007 31.6919 . 0.0000 0.0000 
0.800 0.0173 -0.8596 -5.2198 0.0000 . 0.0000 
0.900 -0.0377 -0.1419 -17.0744 0.0000 . 0.0000 
1.000 -0.0191 0.3837 -2.2115 0.0000 0.0000 
1.100 0.0132 0.1774 7.6918 0.0000 0.0000 
1.200 0.9131 -0.1393 3.2721 0.0000 0.0000 
1.300 -0.0028 -0.1264 =2:9219 0.9000 0.0000 
1.400 -0.0073 0.0332 -2.4304 0.0000 0.0000 
1.500 -0.0006 0.0716 0.7586 0.0000 0.0000 
1.600 0,00534 0.0037. - 1.4098 0.0000 0.0000 
1.700 0.0013 -0.0342 0.0212 0.0000 0.0000 
1.800 -0.0013 -0.011k * -0.6885 0.0000 0.0000 
1.900 -0.0010 0.0137 -0.2033 0.0000 0.0000 
2.000 0.0004 0.0094 0.2822 0.0000 0.0000 


Daljnje smanjenje intervala integracije više nema smisla (osim možda u nekim specijalnim 
slučajevima), jer je dobiveno poboljšanje ustvari neznatno. 





sl. 22 


= 100 = 


Za slijedeći primjer (sl. 23) smanjen je interval integracije za još deset puta, tako da 
iznosi 0.0004 (odnosno potrebno je provesti 5 000 integracija da bi se izračunala cijela krivulja). 
Odštampat će se samo svaka dvijestotineipedeseta vrijednost. 


( 0.0004 ) INTEGRATION INTERVAL 


( 2. ) TOTAL TIME 


TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1,0000 0.0000 -400,0000 0.0000 0.0000 
0,100 -0,3181 10.0489 14.8370 0.0000 0.0000 
0.200 0.4037 2.9241 201.4079 0.0000 0.0000 
0.200 0.2703 -4. 0614 49.2875 0.0000 0.0000 
0.400 -0.1231 =2.5B21 -85.4617 0.0000 0.0000 
0.500 -0.1679 1.3319 -48.3216 0.0000 0.0000 
0.600 0.0135 1.6377 28.5608 0.0000 0.0000 
0.700 0.0866 -0.1965 31.8538 0.0000 0.0000 
0.800 0.0177 -0.8628 -5.1311 0.0000 0.0000 
0.900 -0.0377 -0.1456 -17.1490 0,09000 0.0000 
1.000 -0.0193 0.3845 -2.289k C.0000 0.0000 
1.100 0.0132 0.1798 7.8126 0.0000 0.0000 
1.200 0.0132 -0.1391 3.3245 0.0000 0.0000 
1.300 -0.0028 -0.1278 -2.9180 0.0000 0.0000 
1.400 -0.0073 0.0327 -2.4594 0.0000 0.0000 
1.500 -0.0007 0.0722 0.7477 0.0000 0.0000 
1.600 0.003t4 0.0042 1.4229 9.0000 0.0000 
1.700 0.0013 -0.034L 0.0312 0.0000 0.9000 
1.800 -0.0013 -0.0117 -0.6928 0.0000 0.0000 
1.900 -0.0010 0.0137 -0.21#1 0.0000 0.0000 
2.000 0.0004 0.0096 0.2826 0.0000 0.0000 


Razlika u numeričkim vrijednostima prema prethodnom slučaju postoji samo u četvrtoj deci- 
mali (tj. za cca 0.2%), dok se poboljšanje u obliku krivulje ne može više niti zamijetiti. 


Medjutim, osjetna razlika 
nastala je u vremenu potrebnom za 
proračun. Dok je za prethodne slu- 
čajeve bilo dovoljno cca 2 do 2,5 
minute za računanje i ispisivanje 
rezultata, za ovaj posljednji slu- 
čaj bilo je potrebno oko 9 minuta. 


Iz ovog razmatranja slije- 

di da će numerički rezultati i oblik 
krivulje biti posve zadovoljavajući, 
ako se za interval integracije uzme 
vrijednost koja je 200 do 400 puta 
manja od ukupnog vremena. Pritom se 
ne smije zanemariti niti oblik kri- 
vulje, koji može takodjer diktirati 
izbor intervala integracije. 








-1 
0 


0.2 
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: Kakav će utjecaj na oblik 
krivulje, kao i na rezultat imati 
loše odabrano ukupno vrijeme? Na 
slici 24 prikazana je krivulja koja 
je nacrtana uz premalo, a na slici 
25 krivulja uz preveliko ukupno vri- 
jeme. Jasno je da će se uz premalo 
ukupno vrijeme nacrtati (i izračuna- 
ti) samo dio krivulje koji je premali 
za uobičajeno proučavanje prijelazne 
pojave odziva. Uz preveliko pak ukup- 
no vrijeme, prijelazna će pojava za- 
vršiti suviše brzo. 


Na temelju gornjeg izlaga- 
nja izgleda da bi za prvo računanje 
bilo bolje odabrati veće ukupno vri- 
jeme, pa onda kod ponovnog računanja 
upisati iznos ukupnog vremena kojim 
će se prikazati upravo onaj dio kri- 
vulje kojeg se želi promatrati. To 
je u principu u redu, iako u sebi 
krije nekoliko zamki. Jedna je od 
njih ta, da se može dogoditi da oda- 
brani interval integracije bude pre- 
mali, i da s obzirom na krivulju s 
uskim i oštrim ekstremima računalo 
ne nacrta (i na izračuna) krivulju 
dovoljno točno. (Vidi prethodne pri- 
mjere kad je interval integracije 
premali). 


Druga opasnost leži u inter- 
valu štampanja koji je odabran po 
istoj logici kao i u prvom primjeru 
(za sl. 15). Ako bi se promatrale 
samo numeričke vrijednosti izlaza 3, 
onda bi se na izgled dobio odziv 





Sl. 25 





koji bi imao samo jedno nadvišenje 

s maksimalnim iznosom 0.0011, što 
uopće ne odgovara stvarnosti. No, za 
izlaz 3 nacrtana je i krivulja, pa 
se iz nje može vidjeti kako krivulja 
izgleda i koje su samo vrijednosti 
odštampane (označene križićem). 


Kad bi se medjutim, na te- 
melju numeričkih podataka za izlaz 2 
htjelo procijeniti koje podatke tre- 
ba dati računalu da se krivulja na- 
crta barem približno u redu, onda bi 
se čovjek jako prevario. Iz vrijed- 
nosti izlaza 2 vidi se da se krivu- 
lja nalazi izmedju -0.0005 i +1.5. 
Ukoliko se kao vremenska os želi os- 
taviti srednja crta u okviru crteža, 
bilo bi oportuno postaviti za mini- 
malnu i maksimalnu vrijednost za 
Y-os iznose -2, i +2. . Pogledaju li 
se numerički podaci za izlaz 2 bilo 
kojeg prethodnog proračuna, onda se 
vidi da u vremenu od 0. do 0.5 s 
krivulja ima pozitivni maksimum veći: 
od 10. i negativni veći od -4.5, pa 
su prema tome i podaci za Y-os vrlo 
loše odabrani. 


Ipak, ne može se uvijek re- 
ći da je npr. za sliku 24 ukupno 
vrijeme loše odabrano. Ukoliko se 
želi istražiti samo oblik uzlaznog 
dijela krivulje i prvog nađvišenja, 
onda je ukupno vrijeme od 0.4 s 
upravo odlično odabrano. 
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Sl. 26 


Kako se nakon proračuna i 
crtanja slike 23 može korektno dobi- 
ti samo prvi dio krivulje (tj. kri- 
vulja na slici 24)? 


Kao prvo, treba uključiti 
preklopke 4, 5, 7,8, ilo. 


Iz prethodnih krivulja vid- 
ljivo je da prvo nadvišenje završava 
prije 0.4 s. Dakle, ukupno vrijeme 
neka bude 0.4 . 


U skladu s ukupnom vremenom 
treba odabrati i interval integraci- 
je. Budući da je u tom dijelu krivu- 
lja "vrlo pristojna", interval inte- 
gracije 400 puta manji od ukupnog 
vremena, tj. 

0. 


200 7 0.001 
si 


m 


sasvim će gurno dati dobar rezul- 


tat. 


Buđući da je promijenjeno 
ukupno vrijeme, potrebno je novu 
vrijednost upisati u podatke za X- 
-os, ako se želi da se izračunata 
krivulja rastegne na cijeli format 
okvira. Prema tome, kad računalo 
zatraži podatke za X-os (a zatražit 
će ih zato jer je bila uključena 
preklopka 8), u zagradu za BLOCK 
FOR X-AXIS treba upisati broj 76, u 
zagrađu za MINIMUM VALUE broj 0., i 
u zagradu za MAXIMUM VALUE broj 0.4. 
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SWITCHES SET ON WERE kh 57 8 
( 0.001 ) INTEGRATION INTERVAL 


( ) TOTAL TIME 
0.4 


( ) BLOCK FOR X-AXIS ( ) MINIMUM VALUE ( ) MAXIMUM VALUE 
76 0. 0.4 


(0.02 ) PRINT INTERVAL : 
TIME OUTPUT 3 OUTPUT 2 OUTPUT _ 1 OUTPUT 0 OUTPUT 0 


0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.020 -0.9621 3.6806 -332.6819 0.0000 0.0000 
0.040 -0.8578 6.6031 -249.7657 0.0000 0.0000 
0.060 -0.7038 8.6508 -159. 0869 0.0000 0.0000 
0.080 -0.5179 9.7874 -68.7557 0.0000 0.0000 
0.100 -0.3181 10.0488 14.8456 0.0000 0.0000 
0.120 -0.1211 9.4309 86.3219 0.0000 0.0000 
0.140 0.0589 : 8.3751 141.9752 0.9000 0.0000 
0.160 0.2108 6.7510 179.7677 0.0000 0.0000 
0.180 0.3270 4.9407 199.254 0 0.0000 0.0000 
0.200 0.4037 2.8236 201.4040 0.0000 0.0000 
0.220 0.4403 0.8633 188.3445 0.0000 0.0000 
0.240 0.4395 -0.9025 163. 0491 0.0000 0.0000 
0.260 0.4062 -2.3686 129. 0059 0.0000 0.0000 
0.280 0.3472 -3.4658 89.8915 0.0000 0.0000 
0.300 0.2703 -4.161L4 49.2733 0.0000 0.0000 
0.320 0.1835 -k.4569 10.3600 0.0000 0.0000 
0.340 0.0945 -4.3832 -24.1895 0.0000 0.0000 
0.360 0.0102 -3.9943 -52.3902 0.0000 0.0000 
0.380 -0.0637 -3.3607 -72.9895 0.0000 0.0000 


0.400 -0.1231 -2.5616 -85.4631 0.0000 0.0000 





SWITCHES SET ON WERE kB 7 8 
( 0.025 ) INTEGRATION INTERVAL 


( 10. ) TOTAL TIME 


(76) BLOCK. FOR X-AXIS ( 0. ) MINIMUM VALUE ( 10. ) MAXIMUM VALUE 


( 0.5 ) PRINT INTERVAL 


TIME OUTPUT _ 3 OUTPUT 2 OUTPUT _1 OUTPUT _ 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.500 -0. 1474 1.4938 -37.8331 0.0000 0.0000 
1.000 -0.0106 0.3615 0.8847 0.0000 0.0000 
1.500 0.0011 0.0500 1.1636 0.0000 0.0000 
2.000 0.0005 0.0030 0.2594 0.0000 0.0000 
2.500 0.0001 -0.0005 0.0336 0.0000 0.0000 
3.000 . 0.0000 -0.0002 0.0016 0.0000 0.0000 
3.500 0.0000 -0.0000 -0.0005 0.0000 0.0000 
4.000 -0.0000 -0.0000 -0.0001 0.0000 0.0000 
4,500 -0.0000 0.0000: -0.0000 0.0000 0.0000 
5.000 -0.0000 0.0000 -0.0000 0.0000 0.0000 
5.500 -0.0000 0.0000 0.9000 0.0000 0.9000 
6.000 0.0000 0.0000 0.0000 0.0000 0.0000 


RUN TERMINATED BY SHITCH 0 


AFTER SELECTING DESIRED OPTION PRESS START 


05 


05 
0.2 





5.6 
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Budući da je bila uključena 
preklopka 7, nacrtat će se novi ok- 
vir za sliku, a nakon toga zbog uk- 
ljučene preklopke 5, treba upisati 
interval štampanja. Ako se želi do- 
biti numeričke vrijednosti iz 20 to- 
čaka, treba upisati 0.02 prema 


ukupno vrijeme _ i = 0.02 


20 2 


Sve ostale podatke računalo 
je zapamtilo od prethodnog primjera 
i nije ih potrebno ponovno upisivati. 


Sve što je računalo ispisa- 
lo u vezi ovog proračuna (i slike 24) 
nalazi se na strani 105. 


Za slijedeći primjer (slika 
25 i proračun na strani 106) kod ko- 
jeg je uzeto ukupno vrijeme pet puta 
veće nego u prvom proračunu, pre- 
klopke su mogle ostati uključene još 
cd prethodnog proračuna, budući da 
će se zatražiti promjene istih veli- 
čina. 

400 puta manja vrijednost 
od ukupnog vremena (10. s) dat će 
interval integracije (0.025), a 20 
puta manja vrijednost, interval štam- 
panja (0.5). U podatke za X-os treba 
za broj bloka upisati 76 (vrijeme 
kao varijabla), te 0. za minimalnu i 
ukupno veijeme (10.) za maksimalnu 
vrijednost. 


=. 08. = 


Iz oblika krivulje i iz numeričkih podataka vidjelo se je za vrijeme računanja, da su sve 
vrijednosti svih veličina pale na nulu nakon 4 sekunde odziva. Budući da nije imalo smisla ispisi- 
vati nule do kraja ukupnog vremena, proračun je prekinut uključivanjem preklopke 0. 


Pripremajući se za siijedeći primjer, mogu se odredjene preklopke kojima se mogu zatraži- 
ti pojedine promjene, uključiti već u fazi računanja prethodnog primjera, da bi se uštedjelo na 
vremenu. Tada se može odmah nakon faze računanja pritisnuti tipka PROGRAM START i započeti s una- 
šanjem novih podataka za slijedeći proračun. Računalo uvijek najprije ispiše koje su preklopke bi- 
le uključene, npr.: A 


SWITCHES SET ON WERE 457 8 39 


Medjutim, umjesto da nakon toga zatraži ispisivanje novih podataka ispisivanjem odgovarajućih za- 
grada, može se dogoditi da ispiše 


SWITCHES SET ON WERE hk 57 8 39 
AFTER SELECTING DESIRED OPTION PRESS START 


U tom slučaju je najvjerojatnije preklopka 0 ostala uključena. Da bi se moglo krenuti s 
programom dalje, treba preklopku 0 isključiti i ponovno pritisnuti tipku PROGRAM START. 


U slijedećem primjeru (slika 26 i proračun na strani 109) pokazano je da se položaj nule 
na osi Y može postaviti na bilo kojoj visini. 


Proučavajući prethodne odzive ovog sustava 2. reda, može netko doći na ideju da za mini- 
malnu i maksimalnu vrijednost za os X stavi na primjer, iznose -1. o 0.2 po, recimo, logici da 
krivulja ide od -1., ali da vjerojatno neće imati maksimalno nadvišenje veće od 20 %. Kako će iz- 
gledati krivulja nacrtana s ovakvim podacima prikazuje slika 26. 


Budući da se želi nacrtati krivulja u vremenu 2 sekunde, treba pomoću preklopke 4 zatra- 
žiti promjenu intervala integracije (upisano: 0.005) i ukupnog vremena (upisano: 2.), pomoću pre- 
klopke 5 promjenu intervala štampanja (upisano: 0.1), pomoću preklopke 7 crtanje novog okvira za 
sliku, pomoću preklopke 8 promjenu podataka za X-os (upisano: blok 76, minimalna vrijednost 0. , 


= 109 - 


(0.005 ) INTEGRATION INTERVAL 


( 2. ) TOTAL TIME 
( 3) BLOCK FOR Y-AXIS ( =L. ) MINIMUM VALUE ( 0.2 ) MAXIMUM VALUE 
(76) BLOCK FOR X-AX!S ( 0. ) MINIMUM VALUE ( 2, ) MAXIMUM VALUE 
( .1 ) PRINT INTERVAL 
TIME OUTPUT _ 3 OUTPUT 2 OUTPUT 1 OUTPUT 0. OUTPUT 0 
0.000 -1.0000 0.0009 -4 60. 0000 0.9000 0.9000 
0.100 -0.3173 10.0432 15.0878 0.9009 9.0000 
0.209 0.4036 2.8079 2C1.1606 0.6000 0,0000 
0.300 0.2691 4.1597 #8.8132 6.0000 9.0000 
0.400 -0.12535 -2.5455 -85.3951 $.0000 0.0000 
0.500 -0.1670 1.3364 -47.5087 9.09C0 9.0000 
0.600 0.0141 1.6269 28.6502 0.0909 3.0000 
5.700 2.0862 -0.2031 31.4983 0.0009 9.0900 
9.800 0.0171 . -0.8578 -5.2716 0.0900 9.C000 
6.900 -0.0376 -0.1398 -17.0315 6.0000 0.0999 
1.000 . -0.0190 0.3833 -2.16653 0.0000 0.C600 
1.100 0.0132 0.1760 7.7799 0.0909 9. 0000 
1.200 0.0130 -0.1395 3.2417 6.0009 9.0090 
1.3500 -0.0029 -0.1257. > -2.9241 0.5009 0.0000 
1.400 -0.0072 0.0335 -2.4136 0.0090 0, 0005 
1.500 -0.0006 0.0712 0.7649 0.0000 0.0000 
1.600 0.0034 0.0034 1.4022 9.0000 0.0963 
1.700 9.0012 -0.0341 0.0154 6.0000 0,00050 
3.809 -0.0013 -0.0112 -0.6861 0.0000 9.400) 
1.990 -9.9010 0.0137 -0.1994 0.0000 0.0009 


2.599 2.0004 0.0093 0.2820 0.0009 95,9900 





-0.95 


-0975 


1 


o 


0.01 


sl. 28 


110 


maksimalna vrijednost 2.) i pomoću 
preklopke 9 promjenu podataka za Y- 
-os (upisano: blok 3, minimalna vri- 
jednost -1., maksimalna vrijednost 
0.2). Preklopku 10 korisno je (zbog 
uštede vremena) držati stalno uklju- 
čenu. A 


Kao prvo, iz slike 26 vidi 
se da je X-os pomaknuta prema gore 
i da je mjerilo na Y-osi takvo da je 
ukupna Y-os dugačka 1.2 . Takvo mje- 
rilo nije jako prikladno jer zahti- 
jeva nepotrebna preračunavanja ako 
se vrijednosti očitavaju iz krivu- 
lje. Zato treba nastojati, kad god 
je to moguće, da se jedna od posto- 
jećih triju horizontalnih crta pro- 
glasi za X-os, a jedna od postoje- 
ćih vertikalnih crta okvira za Y-os, 
dok je za mjerilo prikladno uzeti 
takve veličine za minimalnu i maksi- 
malnu vrijednost (za X-os i za Y-os) 
koje su djeljive s 5 ili s 10. Na 
taj način olakšano je snalaženje na 
slici, budući da računalo svaki puta 
kod crtanja novog okvira načini na 
jednoj horizontalnoj i jednoj verti- 
kalnoj crti podjelu na deset jedna- 
kih dijelova. 


S time u vezi preporučljivo 
je za minimalnu i maksimalnu vrijed- 
nost upisati iznose na jedan od sli- 
jedećih tri načina: 

a) za minimalnu vrijednost 0. , 
a za maksimalnu vrijednost broj koji 
je sigurno veći od maksimalne ampli- 
tude krivulje; : 








b) za minimalnu vrijednost nega- 
tivni broj koji je sigurno manji od 
najmanje amplitude krivulje, a za 
maksimalnu vrijednost 0. ; 


c) za minimalnu i maksimalnu 
vrijednost jednaki broj, s time da 
broj za minimalnu vrijednost ima ne- 
gativni predznak. 


Slučaj a) prikladan je za 
krivulju koja nema negativnih, slu- 
čaj b) za krivulju koja nema pozi- 
tivnih vrijednosti, a slučaj c) za 
krivulju koja ima i pozitivne i ne- 
gativne vrijednosti. 


Na slici 26 ujedno je pro- 
vjerena tvrdnja, da pisaljka prili- 
kom crtanja krivulje neće nacrtati 
krivulju izvan okvira, nego će osta- 
ti na maksimalnom položaju, dok će 
bubanj vući papir u smjeru vremenske 
osi. 


Ukoliko je neki dio krivu- 
lje posebno interesantan, može ga se 
povećano nacrtati na cijeli format 
okvira. Neka je npr. posebno intere- 
santan samo srednji dio krivulje 
(sl. 27). U tom slučaju potrebno je 
pomoću preklopki 8 i 9 zatražiti sa- 
mo promjenu minimalne i maksimalne 
vrijednosti za X-os (i eventualno i 
za Y-os), a sve ostale podatke osta- 
viti kao u prethodnom primjeru. 


SWHITCHES SET ON. WERE 1.8 9 
( 3) BLOCK FOR Y-AXIS ( -0.5 ) MINIMUM VALUE ( 0.5 ) MAXIMUM VALUE 


(76) BLOCK FOR X-AXIS ( 0.2) MINIMUM VALUE ( Be4 ) MAXIMUM VALUE 


Numerički podaci za pojedine izlaze odštampani su potpuno jednako kao u prethodnom pri- 
mjeru (rezultati na strani 109). To znači da je računalo izračunalo cijeli zadatak od početka do 
kraja ukupnog vremena, dok je unošenjem minimalne i maksimalne vrijednosti (koje se nalaze unutar 
intervala ukupnog vremena) za vremensku os (X-os), zakočeno gibanje papira na crtaču krivulja od 
početka do trenutka označenog s minimalnom vrijednošću i od trenutka označenog s maksimalnom. vri- 
jednošću do kraja ukupnog vremena. Za cijelo vrijeme računanja pisaljka se giba u smjeru Y-osi i 
crta križiće u trenutku ispisivanja numeričkih podataka. 


Još je jedan dio krivulje posebno interesantan. To je dio samog početka krivulje. Četrde- 
set puta povećani dio početka krivulje prikazuje slika 28. (Vremenska os rastegnuta je 100 puta.) 


Obzirom na novo ukupno vrijeme (0.05), smanjen je interval integracije na 0.0001 i interval Štam- 
panja na 0.0025. 


SWITCHES SET ON WERE kO S To a€e 9 


( 0.00005 ) INTEGRATION INTERVAL 


( ) TOTAL TIME 
0.02 
( ) BLOCK FOR Y-AXIS ( ) MINIMUM VALUE ( ) MAXIMUM VALUE 
3 =1« -0.95 
( ) BLOCK FOR X-AXIS Go ) MINIMUM VALUE ( ) MAXIMUM VALUE 


76 0. 0.02 


( ) PRINT INTERVAL 
0.001 M : 

TIME: OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.001 -0.9999 0.1993 —397.1417 0.0000 0.0000 
0.002 *049996 0.3971 -394.224u 0.0000 0:0000 
0.003 -0.9991 0.5935 -391.2490 0.0000 0.0000 
0.004 -0.9984 0.7884 -388.2167 0.0000 0.0000 
0.005 =0,2975 0.9817 -385. 1284 0.0000 0.0000 
0.006 -0.9964 1,41735 -381.9850 0.0000 0.0000 
0.007 -0.9952 1.3637 -378.7874 0.0000 0.0000 
0.008 =0,9937 145523 517545312 0.0000 0.0000 
0.009 -0.9921 1.7392 e312«2351 0.0000 0.0000 
0.010 -0.9902 1.9245 -368.8819 0.0000 0.0000 
0.011 -0.9882 2.1081 +365.4790 > 0.0000 0.0000 
0.012 -0.9860 2.2899 -362.0272 0.0000 0.0000 
0.013 =0.9836 2.4701 =3559.52175 * 0.0000 0.0000 
0.014 -0.9811 2.6485 -354.9811 0.0000 0.0000 
0.015 *949783 2.8250 . *551.3889 0.0000 0.0000 
0.016 -0.9754 2.9998 -347.7520 0.0000 0.0000 
0.017 -0.9723 3.1728 -344,0714 0.0000 0.0000. 
0.018 -0.9691 3.3439 -340.3481 0.9000 0.0000 
0.019 -0.9656 3.5131 -336.5833 0.0000 0.0000 
0.020 -0.9620 3.6804 -332..1780 0.0000 0.0000 


AFTER SELECTING DESIRED OPTION PRESS START 


Ako nisu potrebni numerički podaci iz pojedinih izlaza, i ako se na kvalitetu proračuna- 
vanog sustava može zaključiti iz oblika krivulje, onda se može vrijeme potrebno za simulaciju na 
digitalnom računalu bitno smanjiti. Da bi se to postiglo, treba postaviti da je interval štampanja 
jednak ukupnom vremenu. U tom slučaju odštampat će se samo vrijednosti izlaza na početku i na kra- 
ju računanja, a krivulja će biti nacrtana bez križića (sl. 29). Crtanje krivulje traje koliko i 
samo računanje. U ovom slučaju crtanje krivulje trajalo je cca 35 sekundi. 





= 14 = 


( 0.005 ) INTEGRATION INTERVAL 


( . ) TOTAL TIME 


( ) BLOCK FOR Y-AXIS ( :) MINIMUM VALUE ( ) MAXIMUM VALUE 
3 =1+ 1 


(76) BLOCK FOR X-AXIS Cc 0. ) MINIMUM VALUE ( Žig ) MAXIMUM VALUE 


( 2. ) PRINT INTERVAL 


TIME OUTPUT _ 3 OUTPUT 20 OUTPUT _1 OUTPUT .0 OUTPUT 0 


0.000 -1.0000 0.0000 *-400.0000 0.0000 0.0000 


2.000 0.0004 0.0000 . 0.2820 0.0000 0.0000 


Dosadašnja razmatranja odnosila su se uglavnom na analizu odabiranja raznih vremenskih 
veličina i mjerila za proračun i crtanje krivulja. 


Prilikom eksperimentiranja s postavljenim zadatkom često je puta potrebno obaviti cijeli 
niz promjena parametara, pa čak i konfiguracije sustava da bi se dobili željeni rezultati. U sli- 
jedećih nekoliko primjera prikazan je postupak kod traženja promjene jednog od parametara uz krat- 
ku analizu dobivenih. rezultata. 


Zadatak se može postaviti na slijedeći način: Ispitati odziv sustava (tj. vraćanje mase 
u ravnotežno stanje) mijenjanjem koeficijenta prigušenja amortizera. 


Promjenu podataka o parametrima i početnim uvjetima zatražit će se uključivanjem preklop- 
ke 2. Da bi se bolje vidio oblik pojedine krivulje i da bi se mogla npr. mjerenjem točaka gdje 
krivulja sijeće vremensku os, provesti usporedba odziva za razne faktore prigušenja, svaka je kri- 
vulja nacrtana na posebnu sliku uz koju su dani i pripadni numerički podaci. 


- 44115 - 


Kako će izgledati odziv ovog sustava ako se postavi vrijednost koeficijenta prigušenja 
amortizera D = 5.656, a svi ostali parametri ostanu isti? 


SWITCHES SET OH WERE 2. 7 
INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1 PAR2 PAR3 
KOEF. 6.1) ( 5.656 ) ( 400. ) ( ) 

TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.2265 12.1731 -21.7466 0.0000 0.0000 
0.200 0.6896 3.4189 295.0843 0.0000 0.0000 
0.300 0.3643 -8.2089 99.2858 0.0000 0.0000 
0.400 -0.4171 -4.8808 -194.4564 0.0000 0.0000 
0.500 -0.3915 4.8124 -129.4002 0.0000 0.0000 
0.600 0.2042 5.0280 110.1272 0.0000 0.0000 
0.700 0.3523 -2.2127 128.3985 0.0000 . 0.0000 
0.800 -0.0549 -4.4087 -46,8891 0.0000 0.0000 
0.900 -0.2808 0.4284 -109.8829 0.0000 0.0000 
1.000 2. -0.0375 3.kL411 4.4566 0.0000 0.0000 
1.100 0.2009 0.6438 84.0188 0.0000 0.0000 
1.200 0.0847 -2.4111 20.2407 0.0000 0.0000 
1.200 -0.1276 -1.1621 -57.6008 0.0000 0.0000 
1.400 -0.0996 1.4897 -31.4237 0.0000 0.0000 
1.500 0.0681 1.2937 34.5610 0.0000 0.0000 
1.600 0.0942 -0.7588 33.3761 0.0000 0.0000 
1.700 -0.0249 -1.1876 -16.6584 0.0000 0.0000 
1.800 -0.0779 0.2380 -29.8186 0.0000 0.0000 
1.900 -0.0032 0.9614 4.1728 0.0000 0.0000 
2.000 0.0578 0.0908 23.6322 0.0000 0.0000 


AFTER SELECTING DESIRED OPTION PRESS START 








Iz krivulje na slici 30 vi- 
di se da je sustav postao oscilator- 
niji što je bilo i za očekivati, jer 
je koeficijent prigušenja smanjen. 
Usput se može reći da je iznos koe- 
ficijenta D odabran tako da ispadne 
krivulja odziva s prigušenjem £ = 
= 0.1. : 


Faktor prigušenja E može se 
dobiti i pomoću izraza 





A 
1 L 
Sri in ak, 
č = 
A 
1 1 
1+ sr ln &, 


gdje su Aj, i A2 iznosi dvaju uzastop- 
nih nadvišenja, i podataka dobivenih 
mjerenjem dvaju uzastopnih nadviše- 
nja iz dobivene krivulje. 


Da bi se dobila krivulja za 
O (sl. 31), treba postaviti 


0. Kao rezultat dobit će se kri- 
vulja neprigušenih oscilacija s vla- 
stitom frekvencijom 

u=/d=/200 = 14,14 
i , : 
a«la= >= 
T= 2,25 0,445 s 
Numerički rezultati ovog 
proračuna nalaze se na slijedećoj 
Ii strani. 


sl. 31 


UBI 


Om 


MH 
LJ 
Nje 
lt 
n 
N 
U 





= 18 = 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1 PAR2 PAR3 
( ) ( ) ( ) ( ) 
1 0. 400. 

TIME OUTPUT 3 OUTPUT _2 OUTPUT 1 OUTPUT _ 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.1548 13.9726 -61.9163 0.0000 0.0000 
0.200 ; 0.9522 4.3256 380.8807 0.0000 \ 0.0000 
0.300 0.4496 -12.6351 179.8376 0.9000 0.0000 
0.400 -0.8131 -8.2378 -325.2526 0.0000 0.0000 
0.500 -0.7014 10.0864 -280.5518 0.0000 0.0000 
0.600 0.5961 11.3614 238.4386 0.0000 - 0.0000 
0.700 0.8860 -6.470h 354.4024 0.0000 0.0000 
0.800 -0.3219 -13.3968 -128.7513 0.0000 0.0000 
0.900 -0.9858 2.4238 -394.3045 0.0000 0.0000 
1.000 0.0167 14. 0488 6.6974 0.0000 0.0000 
1.100 Ds911 1.9562 396.4259 0.0000 0.0000 
1.200 0.2901 -13.5449 116.0282 0.0000 0.0000 
1.300 -0.9014 -6.1497 -360.5541 0.0000 0.0000 
1.400 -0.5692 11.6427 -227.6633 0.0000 0.9000 
1.500 D.7255 9.7548 290.1177 0.0000 0.0000 
1.600 0.7938 -8.6243 > 317.5062 0.0000 0.0000 
1.700 . -0.4796 -12.4259 -191.8591 0.0000 0.0000 
1.800 -0.9424 4.7785 -376.9412 0.0000 0.0000 
1.900 0.1880 13.9068 75. 1884 0.0000 0.0000 


2.000 1.0007 -0.4738 "400. 264L 0.0000 0.0000 





Na 
Na 


.Sl 32 


Za krivulju odziva (sl. 32) 
dobivenu uz £& = 0,5 ( uz uvrštenu 
vrijednost D = 28.28) može se reći 
da praktički ima samo jedno nadviše- 
nje s aperiodskim približavanjem ko- 
načnoj vrijednosti s donje strane. 


Iz numeričkih podataka na 
strani 120 vidi se da postoji i dru- 
go nadvišenje iznosa cca 0,4 %, ali 
ovdje već dolazi u pitanje i tečnost 
računanja (vidi analizu utjecaja 
vremena integracije na točnost re- 
zultata). 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1I PAR2 PAR3 
C.J ( ) ( ) ( ) 
1 28.28 400. 

TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0,4345 7.5661 40.1827 0.0000 0.0000 
0.200. 0.0975 2452171 110.4536 0.0000 0.0000 
0.300 0.1379 -0.9913 27. 1446 0.0000 0.0000 
0.400 0.0224 -0. 9441 -17.7274 0.0000 0.0000 
0.500 -0.0260 -0.0749 -12.5058 0.0000 0.0000 
0.600 -0.0141 0.2040 0.1233 0.0000 0.0000 
0.700 0.0016 0.0863 3.0749 0.0000 0.0000 
0.800 0.0040 -0.0207 0.9971 0.0000 0,0000 
0,900 0.0009 -0.0278 -0.4129 0.0000 0.0000 
1.000 -0.0006 -0.0043 . -0.3798 0.0000 0.0000 
1.100 -0.0004 0.0053 -0.0267 0.0000 0.0000 
1.200 0.0000 0.0028 0.0832 0.0000 ' 0.0000 
1.300 0.0001 -0.0004 0.0343 0.0000 0.0000 
1.400 0.0000 -0.0008 -0.0087 0.0000 0.0000 


' 1.500 -0.0000 -0.0002 -0.0112 0.0000 0.0000 
RUN TERMINATED BY SWITCH 0 








LL 
LILI 


Sl. 33 


: Isto tako za & = 0,707 kri- 
vulja odziva (dobivena za D = 40.), 
trebala bi imati samo jedno nadviše- 
nje s aperiodskim približavanjem ko- 
načnoj vrijednosti s gornje strane 
(sl. 33). Može se reći da je to i 
ispunjeno uz točnost veću od 0,2 % 
uz primijenjeni interval integracije. 


Numerički rezultati koji 
odgovaraju ovom slučaju nakaza se 
na strani 122. 


= 122 = 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK .IC/PAR1 PAR2 PAR3 
( ) ( ) ( ) ( 
ik 40. 400. 

TIME OUTPUT _ 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.5080 6.1823 4.1086 0.0000 0.0000 
0.200 -0.0669 2.4586 71.5784 0.9000 0.0000 
0.300 0.0420 0.1426 22.5073 0.0000 0.0000 
0.400 0.0257 -0.275L4 -0.7188 0.0000 0.0000 
0.500 0.0046 -0.1288 -3.3264 0.0000 0.0000 
0.600 -0.0017 -0.0140 -1.2258 0.0000 0.0000 
0.700 -0.00153 0.0118 -0.0381 0.0000 0.0000 
0,800 -0.0003 0.0066 0.1504 0.0000 0.0000 
0.900 0.0001 0.0010 0.0650 0.0000 0.0000 
1.000 0.0001 -0.0005 0.0055 0.0000 0.0000 
1.100 0.0000 -0.0003 -0.0066 0.0000 0.0000 
1.200 -0.0000 -0.0001 -0.0034 0.0000 0.0000 
1.300 -0.0000 0.0000 -0.0004 0.0000 0.0000 
1.400 -0.0000 0.0000 0.9003 0.0000 0.0000 
1.500 0.0000 : 0.0000 0.0002 0.0000 0.0000 
1.600 0.0000 -0.0000 0.0000 0.0000 0.0000 
1.700 . 0.0000 -0.0000 - -0.0000 0.0000 0.0000 
1.800 0.0000 -0.0000 -0.0000 0.0000 0.0000 
1.900 -0.0000 0.0000 -0.0000 0.0000 0.0000 


2.000 -0.0000 0.0000 . 0.0000 0.0000 0.0000 








NE 
LIL 


sl. 34 


Za č =1 krivulja odziva 
(dobivena uz _ D = 56,56) ima tipični 
oblik aperiodskog odziva s približa- 
vanjem konačnoj vrijednosti s donje 
strane (sl. 34 i rezultati na strani 
124). Vrijeme porasta (vrijeme koje 
je potrebno da krivulja naraste od 
10 %& do 90 & svoje maksimalne vri- 
jednosti) u ovom je slučaju nešto 
smanjeno prema prethodnim primjerima, 
ali je to zato odziv s najmanjim vre- 
menom porasta koji nema nadvišenja. 


= 124 - 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1 PAR2 PAR3 
Co.) ( ) ( ) ( ) 
1 56.56 400. 

TIME OUTPUT _ 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.5867 4.9559 39.9649 0.0000 0.0000 
0.200 -0.2263 2.3638 43,1658 0.0000 0.0000 
0.300 -0.0754 0.8628 18.6438 0.0000 0.0000 
0.400. -0.0233 0.2799 6.5169 0.0000 0.0000 
0.500 -0.0069 0.0851 2.0672 0.0000 0.0000 
0.600 -0.0020 0.0248 0.6200 0.0000 0.0000 
0.700 -0.0005 0.0070 0+1793 0.0000 0.0000 
0.800 -0.0002 0.0020 0.0505 0.0000 0.0000 
0.900 -0.0000 0.0005 0.0140 0.0000 0.0000 
1.000 -0.0000 . 0.0001 0.0038 0.0000 0.0000 
1.100 -0.0000 0.0000 0.0010 0.0000 0.0000 
1.200 -0.0000 . 0.0000 0.0003 0.0000 0.0000 
1.300 -0.0000 0.0000 0.0001 0.0000 0.0000 
1,400 -0.0000 0.0000 0.0000 9.0000 0.0000 
1.500 -0.0000 0.0000 0.0000 0.0000 0.0000 
1.600 -0.0000 0.0000 0.0000 0.0000 0.0000 
1.700 -0.0000 0.0000 0.0000 0.0000 0.0000 
1.800 -0.0000 0.0000 0.0000 0.0000 0.0000 
1.900 -0.0000 0.0000 0.0000 0.0000 0.0000 


2.000 -0.0000 0.0000 0.0000 0.0000 0.0000 





Sl. 35 





š Slijedeća tri proračuna s 
odgovarajućim krivuljama odziva za 
E =1,25 (sl. 35 i numerički rezul- 
tati na strani 126), £=1,5 (sl. 
36 i numerički rezultati na strani 
128) i £=2 (sl. 37 i numerički 
rezultati na strani 130) dobivena 
su uvrštavanjem koeficijenata vis- 
koznog prigušenja D = 70,70, 
D=84,80 i D=113,12. 


= HM6 - 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PARI PAR2 PAR3 
( ( ) ( ) ( ) 
: 70.70 400. 

TIME OUTPUT _ 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.N000 
0.100 -0.6376 4.0868 33.9017 0.000 0.0000 
0.200 * -0.3230 2.2593 30.5299 0.0000 0.0000 
0.300 -0.1598 1.1285 15.8786 0.0000 0.0000 
0.400 -0.0788 0.5573 7.8784 0.0000 0.0000 
0.500 -0.0389 0.2748 3.8874 0.0000 0.0000 
0.600 -0.0192 0.1355 1.9169 0.0000 - C.0000 
0.700 -0.0094 0.0668 0.9451 0.0000 0.0000 
0.800 -0.0047 0.0329 0.4660 0.0000 0.0000 
0.900 -0.0023 0.0162 0.2298 0.0000 0.0000 
1.000 -0.0011 0.0080 D41133. 0.0000 0.0000 
1.000 -0.0006 0.0039 0.0559 0.0000 0.0000 
1.200 -0.0003 0.0019 0.0275 0.0000 0.0000 
1.300 -0.0001 0.0010 0.0136 0.0000 0.0000 
1.400 -0.0001 0.0005 0.0067 9.0000 0,0000 
1.500 -0.0000 . 0.0002 0.0033 0.0000 0.0000 
1.600 -0.0000 0.0001 ' 0.0016 0.0000 0.0000 


RUN TERMINATED BY SWITCH 0 








Prije ovakvog eksperimenti- 
ranja kad se žele mijenjati podaci o 
parametrima, početnim uvjetima, o 
konfiguraciji i generatorima funkci= 
ja, potrebno je provjeriti da li se 
u čitaču kartica nalazi dovoljno 
praznih kartica. Šta znači dovoljan 
broj praznih kartica? To znači da ih 
mora biti barem jedna više nego što 
će biti spomenutih promjena. 


Za promjene u vezi interva- 
la integracije, ukupnog vremena, in- 
tervala štampanja, podataka za X-os 
i Y-os, te za promjene broja eleme- 
nata za čije se izlazne veličine 
traži štampanje numeričkih vrijed- 
nosti, nije potrebno u čitač stav- 
ljati prazne kartice. 
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INITIAL CONDITIONS AND PARANETERS 


IC/PAR NAME BLOCK IC/PARI PAR2 PAR3 
(€...) ( ) ( ) ( 
1- 84.80 400. 

TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.6778 3.5272 27.9998 0.0000 0.0000 
0.200 -0.3972 2.1437 22.9221 0.0000 0.0000 
0.300 -0.2314 1.2508 13.5179 0.0000 0.0000 
0.400 -0.1548 0.7285 1.8773 0.N000 0.0000 
0.500 = -0.0785 0.4243 4.5880 0.0000 0.0000 
0.600 -0.0457 0.2471 2.6721 0.0000 0.9000 
0.700 -0.0266 0.1439 1.5563 0.0000 0.0000 
0.800 -0.0155 0.0838 0.9064 0.0000 0.0000 
0.900 -0.0090 0.0488 0.5279 0.0000 0.0000 
1.000 -0.0053 0.0284 0.3975 0.0000 0.0000 
1.100 -0.0031 0.0166 0.1791 0.0000 0.0000 
1.200 -0.0018 0.0096 0.1043 0.0000 0.0000 
1.300 -0.0010 0.0056 0.0607 0.0000 0.0000 
1.400 -0.0006 0.0033 0.0354 0.0000 0.0000 
1.500 -0.0004 0.0019 0.0206 0.0000 0.0000 
1.600 -0.0002 0.0011 0.0120 0.0000 0.0000 
1.700 -0.0001 0.0006 0.0070 0.0000 0.0000 
1.800 -0.0001 0.0004 0.0041 0.0000 0.0000 
1.900 -0.0000 0.0002 0.0024 0.0000 0.0000 


2.000 -0.0000 0.0001 0.0014 0.0000 0.0900 
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. Treba napomenuti da su 
slike odziva br. 17 - 25 129 - 
38 nacrtane ujedno u tzv. normi- 
ranom obliku, jer se početna i ko- 
načna vrijednost krivulje odziva 
razlikuju za jedinicu. Na taj na- 
čin mogu ove krivulje, kao i nji- 
hovi numerički podaci poslužiti za 
proučavanje krivulja odziva bilo 
kojeg sustava 2. reda. Tako se mo- 
že mjerenjem krivulja odrediti npr. 
iznos pojedinog nadvišenja u pos- 
tocima za pojedine faktore priguše- 
nja, kao i sve ostale karakteristi- 
čne veličine odziva sustava 2, reda, 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1 PAR2 PAR3 
C._2 ( ) C ) ( 
1 113:12 400. 

TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0, 0000 
0.000 0.7371 2.7728 18.8167 0.0000 0.0000 
0.200 -0.5049 1.9134 14.4916 0.0000 0.0000 
0.300 -0.3456 1.3099 9.9292 0.0000 0.0000 
0.400 -0.2366 0.8967 6.7971 0.0000 0.0000 
0,500 -0.0620 0.6139 4.6530 0.0000 0.0000 
0.600 -0.1109 0.4202 3.0853 0.0000 0.0000 
0.700 -0.0759 0.2877 2:0805 0.0000 0.0000 
0,800 -0.0520 0.1969 1.4927 0.0000 0.0000 
0.900 -0.0356 0.1348 1.0218 0.0000 0.0000 
1.000 -0.0243 0.0923 0.6995 0.0000 0.0000 
1.100 -0.0167 0.0632 0.4788 0.0000 0.0000 
1.200 -0.011L4 0.0432 0.3278 0.0000 0.0000 
1.300 -0.0078 0.0296 0.2244 0.0000 0.0000 
1.400 -0.0053 . 0.0203 0.1536 0.0000 0.0000 
1.500 -0.0037 0.0139 0.1052 0.0000 0.0000 
1.600 . -0.0025 0.0095 0.0720 0.0000 0.0000 
1.600 -0.0017 0.0065 0.0493 0.0000 0.0000 
1.800 -0.0012 0.0045 0.0337 0.0000 0.0000 
1.900 -0.0008 0.0030 * (0.0231 0.0000 0.0000 


2,000 -0.0006 0.0021 0.0158 0.0000 0.0000 
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Na slici 38 prikazano je 
sedam krivulja odziva dobivenih mi- 
jenjanjem parametra D pomoću buše- 
nih kartica. (Krivulje su nacrtane 
za slijedeće vrijednosti koeficije- 
nata D: 14.14 0. 5.656 28.28 
40. 56.56 i 113.12). 


Za takav način eksperimen- 
tiranja pripremljen je slijedeći pa- 
ket bušenih kartica: 


- početne kartice 

- kartice s podacima o konfigu- 
raciji 

- jedna prazna kartica 

- kartice s podacima o paramet- 


- dvije prazne kartice 
- kartica s novim parametrom 
- jedna prazna kartica 
- kartica s novim parametrom 
- jedna prazna kartica 


itd. 


Da bi crtanje krivulja išlo 
brže i da se ne štampaju numerički 
podaci za cijeli niz točaka, upisan 
je interval štampanja koji je jednak 
ukupnom vremenu. Od svih izračunatih 
vrijednosti odštampane su samo prva 
i posljednja. 


Nakon što je nacrtana prva 
krivulja, treba uključiti preklopku 
2 (preklopku 2 može se uključiti i 
ranije, ali najranije tek kad raču- 
nalo zatraži podatak o intervalu 
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integracije) i pritisnuti tipku PROGRAM START. Računalo će očitati slijedeću karticu s novim pera- 
metrom. Za vrijeme ispisivanja podataka za novi parametar treba isključiti preklopku 2, jer bi u 
suprotnom, računalo odštampalo četiri para zagrada očekujući da se podaci u vezi još nekih parame- 
tara unesu putem tastature. 


Nakon toga pisaljka se sama postavlja na početno mjesto i počinje crtati novu krivulju, 
Postupak se ponavlja toliko puta koliko se krivulja želi nacrtati. 


Na ovoj i slijedećim stranicama prikazana je cijela lista sa svim podacima u vezi svih 
krivulja nacrtanih na slici 38. 


. CONTINUOUS SYSTEM MODELING PROGRAM 
A DIGITAL ANALOG SIMULATOR PROGRAM FOR THE IBM 1130 


CONFIGURATION SPECIFICATION 


OUTPUT NAME BLOCK TYPE INPUT 1 INPUT 2 INPUT 3 
SILA g W 2 3 0 
BRZINA 2 1 0 1 0 
PUT 3 i U 2 0 0 


INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME , BLOCK IC/PARI PAR2 PAR3 
KOEFICIJENTI 1 14. 1400 400. 0000 0.0000 
MASA 2 0.0000 -0.5000 0.0000 
POCETNI POMAK 3 -1.0000 0.0000 0.0000 
( ) INTEGRATION INTERVAL 

0.005 
( ) TOTAL TIME 
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( ) BLOCK FOR Y-AXIS C ) MINIMUM VALUE ( ) MAXIMUM VALUE 
3 1 


( ) BLOCK FOR X-AXIS ( < ) MINIMUM VALUE ( ) MAXIMUM VALUE 


PREPARE PLOTTER AD PRESS START 
SET PEN ABOUT ONE INCH FROM RIGHT MARGIN 


( ) PRINT INTERVAL 
; TIME OUTPUT( ) OUTPUT( ) QUTPUT( ) OUTPUT( ) OUTPUT( ) 
3 2 Na! 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 


2,000 0.0004 0.0093 0.2820 0.9000 0.0000 
AFTER SELECTING DESIRED OPTION PRESS START ' 


SWITCHES SET ON WERE 2 
INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1 PAR2 PAR3 


M ) < ) 
1 0. 400. 
t + 4“ JN > M: ) 
TIME OUTPUT _ 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1,0000 0.0000 -400.0000 0.0000 0.9000 
2.000 1.0007 -0.4738 400.264 0.0000 0.0000 


AFTER SELECTING DESIRED OPTION PRESS START 
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SWITCHES SET ON WERE 2 
INITIAL CONDITIONS AND PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1 PAR2 PAR3 
Co) ( ) ( ) ( 
1 5,656 400. : 
TIME OUTPUT _ 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT _ 0 
0.000 . -1.0000 0.0000 -400.0000 0.0000 0.0000 
2.000 0.0578 0.0908 23.6322 0.0000 0.0000 


AFTER SELECTING DESIRED OPTION PRESS START 
.SWITCHES SET ON WERE 2 : ' 
INITIAL CONDITIONS AHD PARAMETERS 


IC/PAR NAME BLOCK IC/PAR1I PAR2 : PAR3 
(>) ( ) £ ) ( 
i 28.28 400. 
G “2 ( ) ( ) ( ) 
TIME OUTPUT 3 ' OUTPUT 2 OUTPUT _1 OUTPUT _ 0 OUTPUT _ 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
2.000 *: -0.0000 -0.0000 -0.0003 0.0000 0.0000 


AFTER SELECTING DESIRED OPTION PRESS START 


SWITCHES SET ON WERE 2 
INITIAL CONDITIONS AND PARAMETERS 
IC/PAR NAME BLOCK IC/PARI * PAR2 PAR3 
( 


) < ) < € ) 
1 40. 400. 








TIME 
0.000 
2.000 


AFTER SELECTING 


SVII TCHES SET ON 


IC/PAR NAME 


TIME 
0.000 
2.000 


AFTER SELECTING 


SWITCHES SET ON 


IC/PAR NAME 


TIME 
0.000 
2.000 


AFTER SELECTING 


OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 
-1.0000 0.0000 -400.0000 0.9000 
-0.0000 0.0000 0.0000 0.0000 


DESIRED OPTION PRESS START 


WERE 2 
INITIAL CONDITIONS AND PARAMETERS 


BLOCK IC/PARI PAR2 PAR3 
() ( ) ( ) ( 
I 56.56 400. 
OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 
-1.0000 0.0000 -400.0000 0.0000 
-0.0000 0.0000 0.0000 0.0000 


DESIRED OPTION PRESS START 


WERE 2 
INITIAL CONDITIONS AND PARAMETERS 


BLOCK IC/PAR1I PAR2 PAR3 
Go) ( ) ( ) ( 
1 113.12 400. 
OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 
-1.0000 0.0000 -400.0000 0.0000 
-0.0006 0.0021 0.0158 0.0000 


DESIRED PTION PRESS START 


OUTPUT 0 
0.0000 
0.0000 


OUTPUT _ 0 
0.0000 
0.0000 


OUTPUT 0 
0.0000 
0.0000 





Interesantno je pogledati 
kako izgledaju izlazi iz elementa 2 
(brzina) i elementa 1 (sila). Kri- 
vulje su dobivene uz slijedeće po- 
datke: 


konstanta prigušenja D = 14.14 
interval integracije 0.005 


ukupno vrijeme 2. 
Y-os blok 2 min -15. max 15. 
interval štampanja 2. 


Kako se mijenja brzina ti- 
jekom prijelaznog procesa prikazuje 
slika 39, a promjenu sile prikazuje 
slika 40. , 











Kod slike 39 pada u oči da 


je oblik krivulje brzine ustvari pri- 


gušena sinusoida, dok su krivulje pu- 
ta prigušene kosinusoide. Uz pretpo- 
stavku da je konstanta prigušenja 
jednaka nuli, sinusoida i kosinusoida 
ne bi bile prigušene, pa to vodi na 
ideju da se pomoću dva integratora 
može sastaviti shema koja će uz pri- 
kladno odabrane početne uvjete i pa- 
rametre, generirati funkciju sinusa 
i kosinusa, bez potrebe da se za 
njih sastavlja posebni specijalni 
element. 


Kad se raspolaže s funkci- 
jom sinusa i kosinusa, nije više pro- 
blem dobiti i funkciju tangensa od- 


nosno kotangensa. 


Si. 41 





j Slika 41 prikazuje izlaz, 
tj x(t) kao funkciju brzine v(t) 
tijekom prijelaznog procesa. Kao što 
je i bilo za očekivati dobivena je 
spirala, budući da je izlaz iz in- 
tegratora prigušena kosinusoida, a 
ulaz u integrator prigušena sinuso- 
ida. Ova je krivulja dobivena uz 
iste podatke kao .slika 39, osim što 
je u podatke za X-os upisano: blok 
2, minimalna vrijednost -15. di mak- 
simalna vrijednost 15. 





Na slici 42 prikazana je 
krivulja puta x(t) kao funkcija 
sile F(t) tijekom prijelaznog pro- 
cesa. Za ovaj slučaj upisano je za 
X-os: blok 1, minimalna vrijednost 
-400. i maksimalna vrijednost 400. 


Na ovaj način ujedno je po- 
kazano kako se mogu crtati pored kri- 
vulja y=f(t)_ikrivulje y=f(x), 
pa čak i onda kada su i y i x funk- 
cije vremena. 





-] 
-400 1) 400 


sl. 42 
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Kako primijeniti preklopke Ii4 i 152? 


Primjer neka bude isti kao za sliku 38 s time da nema unaprijed pripremljenih kartica s 
novim parametrima. 


Ovaj način rada kad se kod računanja upotrebljavaju preklopke 14 i 15,posebno je prikla- 
dan kod eksperimentiranja s dinamičkim sustavima kod kojih je teško unaprijed procijeniti ponaša- 
(nje pojedine izlazne veličine. 


Proračun se započne s nekim pretpostavljenim iznosima parametara (za ovaj slučaj neka bude 
D=14.14). Za vrijeme crtanja krivulje (i štampanja numeričkih vrijednosti pojedinih izlaza), re- 
cimo da nismo posve zadovoljni krivuljom koja se crta«i htjeli bismo pokušati promijeniti vrijed- 
nost parametra D, ne bi 1i dobili bolju krivulju. Međiutim, ono što je već nacrtano ne bismo htje- 
li izgubiti i ponovno crtati od početka u slučaju da «idimo da se s novom vrijednošću parametra 
dobiva još lošiji rezultat. Da bi se stanje u trenutkli prekida (uključivanjem preklopke 0) spremi- 
lo u memoriju, potrebno je prije preklopke 0 uključiti preklopku 14. Iz odštampanih podataka vidi 
se da je preklopka 0 uključena za vrijeme štampanja podataka za vrijeme 0.3 . 


SWITCHES SET ON WERE 13 


TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.000 -1.0000 . 0.0009 >400.0000 0.0000 0.0000 
0.100 -0.3173 10,0432 15.0878 0.0000 0.0000 
0.200 0.4036 2.9079 201.1606 0.0000 0.0000 
0.300 0.2691 -4.1597 48.8132 0.0000 0.0000 


RUN TERMINATED BY SWITCH 0 
AFTER SELECTING DESIRED OPTION PRESS START 


Prije nego li se pritisne tipka PROGRAM START, treba uključiti preklopku 2 koja će omo- 
gućiti promjenu parametra D. Odabrana je i upisana nova vrijednost D = 5.656 . Preklopke 2 i 14 
treba isključiti, a zatim pritisnuti tipku EOF. Sada započinje novi proračun iz početka. Nakon 
nekog vremena vidimo da je krivulja lošija od prve i odlučujemo provesti prvi proračun od točke 
prekida do kraja. U tu svrhu potrebno je uključiti preklopke 0 i 15. 
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; Pritiskom na tipku PROGRAM 
START računalo će nastaviti prora- 
čun i crtanje prve krivulje od toč- 
ke prekida do kraja (sl. 43). 


Osim preklopke 15 trebalo 
je uključiti i preklopku 2 da se 
promijeni iznos parametra D natrag 
na prvobitni iznos što u ovom slu- 
čaju nije učinjeno. Iz numeričkih 
podataka na strani 142 vidi se da 
se je proračun nastavio od točke 
prekida, ali s drugim iznosom para- 
metra D (tj. onim koji je zadnji 
unesen u računalo), što ujedno daje 
krivulju s faktorom prigušenja 
E=0.1. 


SWHITCHES SET ON HERE 2 14 


INITIAL CONDITIONS AND PARAMETERS 
IC/PAR NAME BLOCK IC/PARI PAR2 PAR3 


6 3 ( ) ( ) ( ) 
L 5.656 400. : 

TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 QUTPUT_ 0 
0.000 -1.0000 0.0000 -400.0000 0.0000 0.0000 
0.100 -0.2265 12.1731 -21.7466 0.0000 0,0000 
0.200 0.6896 3.4189 295.1843 0.0000 0.0000 
0.300 0.3643 -8.2089 99.2858 0.0000 0.0000 
0.400 -0.4171 -4.9808 -194.456L4 0.0000 0.0000 
0.500 -0.3915 4.8124 -129.4002 0.9000 0.0000 


RUN TERMINATED BY SWITCH 0 


AFTER SELECTING DESIRED OPTION PRESS START 
SWITCHES SET ON WERE 15 


TIME OUTPUT 3 OUTPUT 2 OUTPUT 1 OUTPUT 0 OUTPUT 0 
0.300 0.2691 -4.1597 84,0045 0.0000 0.1000 
0.305 0.2478 -4.3581 74.4527 0.0000 0,0000 
0.400 =0:1922 -3.5017 -96.6012 0.0000 0.N000 
0.500 -0.2567 2.0498 -90.5095 0.0000 0.9000 
0.600 0.0727 3.2414 47.4177 0.0000 0.0000 
0.700 0.2138 -0.7089 81.1923 0.0000 0.0000 
0.800 0.0053 -2.6406 . -12.8285 0.0000 0.0000 
0.900 =-0.1595 -0.2077 -64.9855 | 0.0000 0.0000 
1.000 -0.0488 1.9306 -8.5890 0.0000 0.0000 
1.100 0.1065 0.6987 46.5368 0.0000 0.0000 
1.200 0.0666 -1.2580 19.4398 0.0000 0.0000 
1.300 -0.0615 -0.8796 - -29.5650 0.0000 0.0000 
1.400 -0.0675 0.7005 -23.0216 0.0000 0.0000 
1.500 0.0274 0.8593 15.803 0.0000 0.9000 
1.600 0.0585 -0.2863 21.7791.— 0.0000 0.0000 
1.700 -0.0042 -0.7276 -5.7867 0.0000 0.0000 
1.800 -0.0452 0.0113 -18.0299 0.0000 0.0000 
1.900 -0.0096 0.5513 -0.7052 0.9000 0.0000 


2.000 0.0314 0.1463 13. 3826 0.0000 0.0000 


Prilikom rješavanja zadataka na digitalnom računalu pomoću programa 1130 CSMP, moguće je 
istovremeno rješavati nekoliko odvojenih zadataka, jasno uz uvjet, da su vremena ( interval inte- 
gracije, ukupno vrijeme i interval štampanja) jednaki. Naravno da će se moći dobiti krivulja samo 
za jedan izlaz, ali se numerički podaci mogu dobiti praktički za po jedan izlaz iz pet različitih 
konfiguracija. Pritom treba paziti i na ograničenja u vezi broja upotrebljenih elemenata. 


Isto je tako moguće, pored proračuna glavnog problema, istovremeno obaviti i nekoliko 
drugih proračuna koji mogu i ne moraju biti u vezi glavnog proračuna.Na primjer, kad bi se kod do 
sada razmatranih slučajeva parametar D mijenjao nasumce, bilo bi zgodno kad bi računalo usput iz- 
računalo i vrijednost faktora prigušenja č. : 


To se može bez problema provesti, ako se zna da je (za ovaj slučaj) iznos faktora prigu- 
šenja definiran izrazom 





Neovisno o slici 15 prema kojoj su radjeni svi dosadašnji proračuni, treba sastaviti blo- 
kovsku shemu za gornju jednadžbu (sl. 44). 





SI. 44 
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Budući da petlja nije zatvorena, nije potrebno upotrijebiti nikakav memorijski element 
(npr. integrator). Osim toga računanje faktora £ nije u ovom slučaju ovisno o vremenu (kao što 
je to na primjer, kod sustava s vremenski promjenljivim koeficijentima), pa će se rezultat ispi« 
sati odmah u prvom redu ispisivanja numeričkih vrijednosti pojedinin izlaza. 

Naravno, da i za ovu blokovsku shemu treba pomoću bušenih kartica (ili pomoću tastature) 
unijeti u računalo podatke o konfiguraciji i podatke o parametrima, - 


Nakon završetka rada na jednom primjeru, potrebno je izvaditi kartice iz čitača kartica. 
Da bi se ispraznio čitač kartica, tj. da bi se dvije kartice koje se nalaze unutar čitača dobile 
van, treba pritisnuti tipku MPRO koja se nalazi na čitaču kartica. 


Ako se nakon prethodnog primjera želi početi s nekim drugim primjerom iz početka, potre- 
bno je pritisnuti tipku INT REQ na centralnoj jedinici, a od početnih kartica staviti samo karti- 
ce // JOB i // XEQ CSMP. Iza njih dolaze odmah kartice s podacima o konfiguraciji. 


Da bi se računalo pokrenulo, potrebno je pritisnuti na čitaču kartica tipku START i na 
centralnoj jedinici tipku PROGRAM START. U tom trenutku zasvijetli RUN i ponovno počne prepisiva- 
nje programa iz diska u centralnu memoriju. 
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PRIMJER 2 


Potrebno je proučiti ponašanje regulacijskog sustava tijekom prijelaznog procesa, zadanog 
prema slici 45. 


regulator 





generator 


SI. 45 


Zadatak je formuliran vrlo jednostavno: potrebno je regulirati brzinu motora tako da br- 
zina odmatanja kabela s bubnja bude konstantna i da iznosi 15 m/s (tj. da se svake sekunde odmota 
s bubnja 15 metara kabela). 
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U ovom slučaju gibanje je kružno, pa umjesto mase dolazi inercija, umjesto sile dolazi 
zakretni moment, a umjesto pravocrtnog ubrzanja dolazi kutno ubrzanje. 


Iako na prvi pogled izgleda jednostavan, ovaj bi zadatak bio teško rješiv bez elektronič- 
kog računala, budući da su koeficijenti promatranog sustava vremenski promjenljivi. 


Da bi se linearna brzina odmatanja kabela održala konstantnom, potrebno je povečavati br- 
zinu motora, jer kako se kabel odmata polumjer bubnja se smanjuje. S druge strane, smanjenjem ko- 
ličine kabela na bubnju, smanjuje se inercija bubnja, a s time i potrebni zakretni moment kojeg 


mora razviti motor. 

Da bi se uopće moglo pristupiti proračunu, potrebno je poznavati svojstva pojedinih ele- 
menata ovog regulacijskog kruga, čiju blokovsku shemu prikazuje slika 46, kao i numeričke vrijed- 
nosti svih parametara. 





SI. 46 


Odmatajući se, kabel prolazi izmedju valjaka čija je brzina okretanja proporcionalna br- 
Zini kabela, Masa,odnosno inercija ovih valjaka može se u proračunu zanemariti. Brzina okretanja 
valjaka mjeri se tahometrom koji ima linearnu karakteristiku i čija vremenska konstanta iznosi 
0.5 s. Prijenosna funkcija tahogeneratora može se dakle prikazati izrazom 


(1) 





FalP) = sz 
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Izlaz iz motora odredjen je pojačanjem motora i vremenskom konstantom. Zbog jednostavno- 
sti uzeta je prijenosna funkcija prvog reda definirana omjerom zakretnog momenta prema ulaznom na- 
ponu: 


V, 
i: M (2) 
PAPE Tr 
M TMP +1 
gdje je: Vi = 500 pojačanje motora 
TM =1 vremenska konstanta motora 


Zakon po kojem se smanjuje polumjer bubnja glasi: 


R=R,- 0.0008 a (3) 
gdje je: R trenutni polumjer bubnja 
R, =4 polumjer punog bubnja 
Ru 2 polumjer praznog bubnja 
[07 ukupni zakret bubnja u radijanima 


Zakon po kojem se mijenja inercija bubnja J, glasi: 
J = 18.5 R* - 221 (4) 
Stvarna brzina kabela _ v može se dobiti kao umnožak trenutnog polumjera i kutne brzine 
bubnja: 
v = R6 (5) 
Ubrzanje bubnja dano je omjerom izmedju zakretnog momenta motora Ma i ukupne inercije J: 
M 
s —mn 
J 


a = 


(6) 


Sada se može pristupiti crtanju simulacione blokovske sheme koja se može nacrtati u is- 
tom rasporedu kao blokovska shema na slici 46. 
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Motor zadan prijenosnom funkcijom (2) može se odmah nacrtati. Za tu svrhu dovoljan je je- 


dan integrator prema slici 47. 
e ž 


SI. 47 


Budući da oba ulaza treba pomnožiti s nekim parametrima, treba ući na drugi i treći ulaz 
integratora. Na taj način, za cijeli motor potrošen je samo jedan element. Elementi simulacione 
blokovske sheme mogu se odmah označiti brojevima, pa neka ovo bude element broj 1. 


Regulator se mora sastojati od takvih elemenata da omogući usporedjivanje stvarne i na- 
mještene (željene) veličine, kao i množenje tog signala pogreške s iznosom pojačanja regulatora. 
Budući da u vezi regulatora nisu postavljeni nikakvi posebni zahtjevi, trebalo bi se odlučiti na 
vrst regulatora kao i odabrati odgovarajuće parametre. Najjednostavnije je uzeti idealni propor- 
cionalni regulator. Faktor pojačanja neka mu bude 1.5 . Simulacionu blokovsku shemu moguće je iz- 
vesti pomoću dva elementa: sklopa za pribrajanje konstante i pojačala prema slici 48. Parametar 
elementa 0 neka predstavlja namještenu vrijednost vo =15,a pojačalo pojačanje regulatora 
ML S La o 

R 





Znatno kompliciranije je sastavljanje simulacione sheme za bubanj s kabelom. Ulaz u bubanj 
je zakretni moment motora, a izlaz neka bude linijska brzina kabela v. Kako početi? Ideju može dati 


izraz (6) koji sadrži veličinu Ma (zakretni moment motora). Ako se Ma podijeli s J, dobit će se 4 
(sl. 49). 
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Integrira li se & dobit će se 4 koji će trebati kasnije prema izrazu (5), pomnožen s R 
dati brzinu v (sl. 50). 
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Iz izraza (3) vidi se veza izmedju veličina a iR. a se može dobiti Podeedastijem &, pa 
shema izgleda (sl. 51): 





Sl. 51 


R, je prvi, a -0.0008 drugi parametar integratora broj 4, pa prema tome s a4 treba ući na 
drugi ulaz integratora. 


Sada kad je dobiven R, moguće je pomoću njega dobiti inerciju J koja je potrebna prema 


slici 49. J se može izračunati prema izrazu (4), kojem odgovara slika 52. 
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Stvarna brzina gibanja kabela _ v dobit će se prema izrazu (5) za što je potrebno jedno 
množilo (sl. 53). 





SI. 53 


Na kraju dolazi još mjerni pretvarač (tahogenerator) čija se simulaciona shema (sl. 54) 
može sastaviti prema prijenosnoj funkciji (1). 





SI. 54 
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Ovdje je upotrebljen mali trik. Poznato je da su regulacijski sustavi, sustavi s negativ- 
nom povratnom vezom. Znači, veličina koja ulazi u zbirno mjesto (element broj 12) mora imati nega- 
tivni predznak. To je najjednostavnije izvesti s negativnim pojačanjem tahogeneratora. Budući da 
tahogenerator ima vremensku konstantu, treba integratoru 10 izvesti povratnu vezu i pomnožiti je 


la=bk sus 
S Tam 5 2.0 


Kako pojačalo 11 već množi izlaz iz integratora 10 s -2., onda se može taj isti iznos 
uzeti za povratnu vezu integratora. Ipak treba definirati drugi parametar integratora 10, jer je 
prvi ulaz u integrator već potrošen za ulazak veličine v. (Drugi i treći ulaz u integrator idu 
preko parametara i ukoliko se oni eksplicitno ne definiraju, računalo ih postavlja na nulu, pa bi 
i odgovarajući ulazi ispali jednaki nuli.) 
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10 
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SI. 56 
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Na taj način dobiveni su 
svi dijelovi simulacione blokovske 
sheme koje treba samo medjusobno 
spojiti, što prikazuje slika 55. 


Podatke iz ove blokovske 
sheme treba prepisati u Listu poda- 
taka o konfiguraciji i Listu poda- 
taka o početnim uvjetima i parame- 
trima i prema njima izbušiti karti- 
ce. Kako su ispunjene ove liste i 
koje su vremenske i ostale veličine 
odabrane za ovaj primjer vidi se iz 
priložene kompletne liste koju je 
računalo odštampalo. Slika 56 pri- 
kazuje vremenski dijagram stvarne 
brzine kabela. 
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CONTINUOUS SYSTEM MODELING PROGRAM 
A DIGITAL ANALOG SIMULATOR PROGRAM FOR THE IBM 1130 


CONFIGURATION SPECIFICATION 


OUTPUT NAME BLOCK TYPE INPUT 1. IHPUT 2 INPUT 3 
MOMENT MOTORA 1 0 13 1 
UBRZANJE BURNJA 2 / 1 8 0 
BRZINA BUBNJA 3 I 2 0 n 
POLUMJER m I 0 3 0 

5 X u m 0 

6 X 5 5 0 
' 7 G 6 0 0 
MOMENT INERCIJE 8 0 7 0. 0 
STVARNA BRZINA 9 X 4 3 0 
TAHOGENERATOR 10 I 9 11 o 
MJERENA BRZINA 11 G 10 0 0 
POGRESKA 12 o) 11 0 0 
POJACANJE 13 G 12 0 0 

INITIAL CONDITIONS AND PARAMETERS 

IC/PAR NAME BLOCK IC/PAR1 2 PAR2 PAN3 
MOTOR 1. 0.0000 500. 0000 -1.09000 
BUBRNJ m 4, 0000 -0.09008 0.0000 

7 18.5000 0.0000 0.0n00 
MOMENT INERCIJE 8 -221.0000 0.0000 0.0000 
TAHOGENERATOR 10 0.0000 1.9000 0.0000 
POJACANJE TAKA 11 -2.0000 0.0000 0.0000 
ZELJENA BRZINA 12 15.0000 0.0000 £.nn90 
POJACANJE REG. 13 1.5000 0.9000 0.0000 


(01 ) INTEGRATION INTERVAL 
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(20.  ) TOTAL TINE 
( 9) BLOCK FOR Y-AXIS — ( 0. ) MINIMUM VALUE (20. ) MAXIMUM VALUE 
( ) BLOCK FOR X-AKIS € ) MINIMUM VALUE ( ) MAXIMUM VALUE 


PREPARE PLOTTER AND PRESS START 
SET PEN ABOUT ONE INCH FROM RIGHT MARGIH 


( , ) PRINT INTERVAL 
TIME OUTPUT( 9) OUTPUT( 4) OUTPUT(11) OUTPUT(12) OUTPUT( 1) 
0,000 0.0000 4.0000 0.0000 15.0000 , 0.0000 
1.000 3.6174 3.9997 -1.6632 13.3368 6815.1718 
2.000 10.2764 3.9984 -7.1118 7.8882 7415.5566 
3.000 15.8109 3.9957 ,, -13.1578 1.8422 4753.0642 
4,#00 18.5094 3,9922 ' -17.1868 -2.1868 1403.6984 
5.000 18.5716 3.9885 -18.5070 -3.5070 -083.0552 
6.000 :17.1607 3.9849 -17.8127 =2.48127 -1887.2790 
7.000 15.5076 3.9816 -16.2910 =1.2910 -1616.4550 
8.000 14.28R0 3.9786 -14.9276 0.0724 -810.2971 
9.000 14.014. 3.9757 -14.2018 0.7982 -85.5717 
10.000 14,2029 3.9729 >: -14.1205 0.8795 408.6589 
11.000 . 14.6245 3.9700 -14.4270 0.5730 491.2933 
12.000 15. 0030 3.9670 -14,9224 0.1776 341.5670 
13.000 15.2034 3.9639 -15.1058 -0.1058 124.0540 
14.000 15.2235 3.9609 :-15:2116 -0.2116 -40.1055 
15.000 15.1358 3.9578 -15.1761 -0.1761 -109,267G6 
16.000 15.0245 3.9547 1540771 -0.0771 -97.2179 
17.000 14.9461 3.2517 -14.9838 0.0162 , -45.3243 
18.000 14.9182 3.9487 -14.9321 0.067. A 7.1809 
19.000 14.9298 3.9457 -14.9248 0.0752 38.3381 
20.000 14.9585 3.9426 -14.9451 0.0549 44,7954 


AFTER SELECTING DESIRED OPTIOHM PRESS START 
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Na temelju odštampanih rezultata i oblika krivulje može se procijeniti koliko odabrani 
regulator zadovoljava ili ne. Sigurno je, na primjer, da motor ne bi u početku izdržao cca 150 
puta veće opterećenje od niminalnog, pa bi i to sigurno bio jedan od razloga zbog kojih bi bilo 
potrebno provesti odredjene zahvate u ovom regulacijskom krugu. 


1) 
2) 
3) 
4) 


5) 
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